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Chapter 1 

Heavy Ion Physics 



1.1 Introduction 

Colliding heavy ions in particle accelerators offers a unique opportunity to study 
the strong interaction of matter in the regime of extremely high densities and 
temperatures. It is believed that in such collisions temperatures and densities 
are reached that prevailed in the universe the first few microseconds after the 
Big Bang. 

In the Standard Model of particle physics, the strong interactions between 
the fundamental quark constituents of matter are described by a field theory 
called Quantum Chromo Dynamics (QCD) 1 . In this theory, the quarks carry 
a strong charge, called color, and the strong force is mediated between the 
colored quarks by the exchange of gluons, which are the quanta of the strong 
field. A very important feature of QCD is that the gluons also carry color 
charge, so that they do not only act as mediators but also themselves couple 
to the strong force. It turns out that, as a consequence, the potential increases 
with increasing distance between the color charges. This is in sharp contrast 
with the field theory of Quantum Electro Dynamics (QED) |2], where the force 
between electrically charged particles is mediated by the electrically neutral 
photon. Here the potential vanishes for large distances. 

The behavior of the strong coupling with varying distance, which is related to 
the behavior of the potential as discussed above, has profound phenomenological 
consequences. First, the coupling between the colored quarks becomes weak at 
short distances, a property called asymptotic freedom. Such short distances are 
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probed in hard scattering processes, where the momentum exchange between 
the participating quarks is large. Since the strong coupHng is weak in the 
hard regime, the interaction cross sections can be calculated in a framework 
called perturbative QCD (pQCD). Because they are calculable, hard scattering 
processes form a unique probe of the constituents of matter, while being, at the 
same time, a testing ground for the validity of QCD. In this way, QCD has 
been firmly established as the correct theory of the strong interaction in the 
last four decades by performing a large variety of experiments on deep inelastic 
scattering of electrons and muons on protons and neutrons and by the study 
of electron-positron and proton- (anti)proton collisions at large centre of mass 
energies in storage rings. 

The strong coupling increases with the distance between the quarks, and 
the interaction becomes, in fact, so strong, that in ordinary matter the quarks 
are permanently confined to colorless hadrons. In this regime of large distances 
or, equivalently, small momentum transfers, pQCD breaks down, so that it 
cannot be used to calculate soft scattering cross sections from first principles. 
However, recently much progress has been made in the understanding of the 
non-perturbative domain by so-called lattice QCD calculations, where the QCD 
field equations are numerically solved on a discrete space-time lattice [3]- 

One of the remarkable results of lattice QCD is the prediction that hadronic 
matter at sufficiently high temperatures and densities will undergo a phase 
transition to a state of quasi-free quarks and gluons. This deconfined dense 
state of matter is called a Quark Gluon Plasma (QGP). In Figure 11.11 are 
shown lattice QCD calculations of the energy density (e) divided by the fourth 
power of the temperature (T) [4] . This dimensionless quantity is proportional to 
the effective number of degrees of freedom available in the medium. Below the 
critical temperature (Tc) the medium consists mainly of confined hadrons, while 
above Tc the quarks and gluons become deconfined, causing a rapid increase in 
the number of degrees of freedom. Figure 11.11 shows that the phase transition 
occurs when nuclear matter is heated to a temperature Tc of about 175 MeV, 
corresponding to an energy density of 0.7 GeY/ixs? . 

In the limit of an ideal Stefan-Boltzmann gas, the equation of state (EoS) of 
a QGP is given by 

PsB-^ and e^^ = ^gT\ (1.1) 

where Psb is the pressure, £sb the energy density, g the effective number of 
partonic degrees of freedom, and T is the temperature [5J . The effective number 



1.1. INTRODUCTION 



16 
14 
12 
10 
8 
6 
4 
2 




RHIC 




1. 



+ + 



LHC 



3 flavour 
2 flavour 



avour 



= (173 +/- 15) MeV 
~ 0.7 GeV/fm^ 



T [MeV] 



300 



400 



500 



600 



Figure 1.1: The energy density e scaled by calculated from lattice QCD with (2,3) de- 
generate quark flavors, as well as with two light and one heavy (strange) quark. The arrow 
on the right-hand side shows the Stefan-Boltzmann limit for an ideal quark-gluon gas. The 
critical temperature Tc and the temperatures which can presumably be reached by RHIC and 
LHC are also indicated. Figure taken from 0], 



of partonic degrees of freedom is given by 

7 



■ 5q +5g, 



(1.2) 



where f/q and gg are the degeneracies of, respectively, the quark and gluon states. 
Each quark flavor has a quark/antiquark state, two spin states, and three color 
states, whereas each gluon has two spin states and eight color states. The total 
degeneracy is, therefore, given by 



7 21 
(;=-xnt x2x2x3 + 2x8 = -z-n{ 

8 2 



16, 



(1.3) 



which yields the value g = 37(95/2) for an rii — 2(3) flavor QGP. This is an 
order of magnitude larger than for a hadron gas, where g ~ 3. 

The horizontal arrow in Figure 11.1! indicates the Stephan-Boltzmann limit 
for a QGP with nf = 3 light flavors. The lattice QCD calculation shows that 
e/T^ above Tc remains far below this limit, indicating that a QGP, according 
to these calculations, does not behave as an ideal gas of quarks and gluons. 
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The possible existence of a QGP was conjectured before the advent of lattice 
QCD calculations, and already in the 1980s experiments started to look for 
signatures of this plasma in heavy ion collisions. This initiated the rapidly 
developing field of heavy ion physics, and led to a large series of experiments 
performed at the AGS in Brookhaven, the ISR and the SPS at CERN, and, since 
the year 2000, at the Relativistic Heavy Ion Collider (RHIC) at the Brookhaven 
National Laboratory (BNL) in the USA. 



1.2 Heavy ion collisions 

To describe a particle collision, we denote by pA the 4-momentum of particle 
A moving along the beam {z axis), and by pb the 4-momentum of particle B 
moving in the opposite direction. The Lorentz-invariant measure of the square 
of the center-of-mass energy available in the collision is 

s^ipA+Psf- (1.4) 
The Lorentz-invariant inclusive cross section of the scattering process 

AB CX 



g "^^ , 7"^^^ = —r^^^, (1.5) 



is defined by 

^d^crj AB CX) 

dp"^ Pt dpT dy d(j) ' 

where C is the final state particle being measured and X denotes all other parti- 
cles produced in the coUison [6| . Because of azimuthal symmetry, it is convenient 
to separate longitudinal and transverse momentum components. In Eq. (jl.Sp . 
E and p are the energy and 3-momcntum, pT is the transverse component of the 
momentum, (f) is the azimuthal angle, and y is the rapidity of particle C in the 
center-of-mass frame. The rapidity is a measure of the longitudinal momentum 
component (pl) and is defined by 

.4>4I^)^ (1.8) 



2 \E-pL 

The rapidity variable has the advantage of being additive under Lorentz boosts 
along the z axis. Another commonly used variable is the pseudorapidity {rj) 
defined by 

= -lntan(6l/2), (1.7) 
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which is simply a measure of the polar angle {9) and does not depend on the 
particle mass. This is, therefore, a convenient variable, since it can be calculated 
without knowing the particle identity. In the limit E — \/ — vn? « p ^ m of 
very energetic particles, the pseudorapidity r\ approaches the rapidity y, because 
particle masses can then be neglected. 

Because atomic nuclei are spatially extended objects, a characteristic of 
nucleus-nucleus collisions is the impact parameter (&), which is the transverse 
distance between the centers of the two colliding nuclei, as shown in Figure [L2l 




Figure 1.2: Transverse view of two colliding nuclei defining the reaction plane, the impact 
parameter 6, and the azimuthal angle </< of a produced particle with respect to the reaction plane. 



Other measures of the collision centrality are the number of participants 
(iVpart) and the number of binary collisions (iVcoii). The number of binary 
collisions A'coU is defined as the number of individual inelastic nucleon-nucleon 
collisions that happened during the nucleus-nucleus collision. The number of 
participants A^part is defined as the number of nucleons that suffered at least 
one inelastic collision with another nucleon. The relation between the impact 
parameter h and the number of collisions iVcoii or A^part is calculable in the 
framework of the Glauber model [7] . 

Experimentally, the centrality of a heavy ion collision is estimated from 
a measurement of one or more quantities that vary monotonically with the 
impact parameter. Such quantities are the charged particle multiplicity (iVch), 
the transverse energy (-Et) of all charged particles emitted near midrapidity, or 
the forward energy {Ep) measured close to the beam line. The relation between 
the observables and the impact parameter is established by Monte Carlo event 
generators that model nuclear collisions at relativistic energies [S]. 
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The range of impact parameters can be represented as a fraction of the total 
geometric cross section. It is customary to define centrahty classes as adjacent 
intervals in b that contain a certain percentile of the differential cross section 
da/db. For instance, a 0-5% centrality class contains events with five percent 
of the smallest impact parameters, such that it corresponds to five percent of 
the total geometric cross section. 



1.3 Heavy ion physics at RHIC 

RHIC is a multipurpose colliding beam facility [HI HO] , capable of accelerating 
protons, deuterons, and heavy ions over a broad energy range. At present, RHIC 
has delivered colliding beams of protons, deuterons, copper, and gold ions with 
beam energies of up to 100 GeV per nucleon [TTl [T^]. 

An estimate of the energy density in the created medium is obtained using 
the Bjorken formula 



where tq is the formation time and R is the initial radius of the expanding 
system. Using the value dEx /drj = 503 ± 2 GeV measured in central Au + Au 
collisions [131 and taking R — 1.2A^^^ fm, together with reasonable guess for 
the parameter value tq = 1 fm/c, an initial energy density of about 5 GeV/fm'^ 
is calculated. This is well above the critical energy density of about 1 GeV/fm'^ 
predicted by lattice QCD for a phase transition to the quark-giuon plasma, as 
shown in Figure fTTTI A major part of the physics program at RHIC is, therefore, 
to measure particle production in high energy nuclear collisions with the goal 
to study the properties of the state of matter (presumably a QGP) produced in 
such collisions. 

Particles emitted with large transverse momentum are important probes of 
the medium produced in the collision, because they most likely originate from 
high energetic partons that propagate through and couple to the created medium 
and thus carry information about its properties. A convenient way to observe 
medium-induced modification of particle production is to compare a nucleus- 
nucleus collision {A + B) with an incoherent superposition of the corresponding 
number of individual nucleon- nucleon collisions {N + N) . This is done via the 
nuclear modification factor (Rab), defined as the ratio of the particle yield in 
nucleus-nucleus collisions and the yield in nucleon-nucleon collisions scaled with 



£Bj - 



dEr 1 
dy ctq 7ri?2 ' 



(1.8) 
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the number of binary collisions A^coii, 

^ d'^NAB/dpTdy 
"^^^ (Tab) d^aP+P/dpTdy- ^''"'^ 

Here (Tab) is the nuclear overlap function that is related to the number of 
inelastic nucleon-nucleon collisions in one A + B collision through 

(T^b) x<,^ = (7Veon). (1.10) 

In the absense of medium effects, the nuclear modification factor is unity, while 
Rab < 1 indicates a suppression of particle production in heavy ion collisions, 
compared to an expectation based on an incoherent sum of nucleon-nucleon 
collisions. 

In Figure fTTSl we show the ratio Raa of charged hadron production, as a func- 
tion of pt, measured by the STAR Collaboration in central Au -I- Au collisions 
at -^/s^vTV = 200 GeV [15] (the quantity is the centcr-of-mass energy of an 
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individual nucleon-nucleon collision). It is evident that charged particle produc- 
tion in Au + Au collisions is significantly suppressed, compared to that in p + p 
collisions at the same center-of-mass energy, in particular at large j»t ~ 8 GeV/c, 
where Raa reaches a value of about 0.2. 

Also shown in Figure 11.31 is the nuclear modification factor measured in 
minimum bias (no centrality selection) and central d + Au collisions. This 
measurement is important to distinguish between initial and final state effects. 
Since we can safely assume that in d + Au collisions no hot and dense medium 
is created, the presence of a suppression would indicate initial state effects, such 
as nuclear modification of the parton densities in the gold nucleus. It is seen 
from Figure [T31 that such suppression is absent in d + Au collisions, indicating 
that the suppression observed in Au + Au collisions is a final state effect caused 
by the dense medium created in such collisions. 

A significant enhancement RdA > 1 seen in d + Au collisions in the region 
2 < pt < 7 GeV/c in Figure 11.31 can be explained by the so-called Cronin 
effect [E]. This effect is likely caused by multiple scattering of the projectile 
parton inside the target nucleus, which acts as an additional transverse momen- 
tum kick of the parton, overpopulating the pr > 2 GeV/c region. Since there 
is an indication in Figure 11.31 of a possible suppression in d -I- Au collisions at 
Pt > 8GeV/c, it is interesting to measure the RdA factor at even higher pT- 
This thesis presents such a measurement. 

For peripheral collisions, the number of participant nucleons is small and the 
creation of a dense medium is not expected. This is illustrated in Figure 11.41 
which shows the centrality dependence of Raa for charged hadrons as measured 
by STAR in Au-I- Au collisions. Indeed, the large suppression observed in central 
collisions gradually vanishes with decreasing centrality. This suggests that, 
instead of p+p interactions, peripheral collisions can be used as a reference. This 
is done through the ratio of particle production in central (C) and peripheral 
(P) events: 

_ {N,on)pd^Nc/dpTdy 
"""^ (A^coii)cd2iVp/dpTdy' ^ ' ^ 

The advantage of this measure is that no p + p reference data are needed. The 
disadvantage is that a stronger model dependence is introduced, because the 
uncertainties in (iVcou) are much larger for peripheral collisions. In Figure [L5] is 
shown Rc p for charged hadrons measured in Au + Au collisions by STAR [T^ . 

To provide a useful reference, it is important to measure particle production 
in nucleus -nucleus interactions, as well as in the p + p collisions, under the same 
experimental conditions. For instance, prior to having the first p + p collisions 
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Figure 1.4: Ratio Raa of charged hadron production, as a function of collision centrality, 
measured by the STAR Collaboration in Au + Au collisions at ^sjvjv = 200 GeV. Figure 
taken from |17) . 



delivered by RHIC, both STAR and PHENIX collaborations have published 
the measurements of Raa [Ml ISj based on p + p and p + p reference spectra 
obtained from a large body of world data, extrapolated to RHIC energies. 
These extrapolations yielded significant systematic uncertainties and more 
precise measurements of Raa [13 HO] only became available when p + p 
reference data were taken at RHIC in 2001-2002. 
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Figure 1.5: Ratio Rcp of the charged hadron production measured by the STAR Collaboration 
in Au + Au collisions at y'sjvjv = 200 GeV. Figure taken from | 17| . 



A detailed study of the intermediate- and high-pT production of various 
hadron species shows that there is a systematic difference between meson and 
baryon production in Au + Au coUisions, as illustrated in Figure flTBl [21] . The 
Rcp ratio for identified hadrons is shown separately for mesons (a) and baryons 
(b) , and the clear difference between them suggests that the particle production 
in this pt range depends not on the mass of the hadron but rather on the number 
of valence quarks contained within it. This can be explained naturally in the 
quark recombination model for hadron formation, rather than fragmentation. 
We do not discuss this model here and refer to [H [H [Ml ISS HI HZ] for 
details. The measurement of Rcp for neutral pions and eta mesons would also 
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(a) and baryons (b) in Au + Au collisions at y'sjvjv = 200 GeV. Figure taken 



be interesting in context of this observation. 

This thesis presents a basehne measurement with the STAR detector of 
neutral pion and eta meson production in p+p and d+Au coUisions at a center-of- 
mass energy of y^s^N = 200 GeV. The neutral pion spectrum complements that 
of the charged pions measured in STAR in the range 0.35 < pt < 10 GeV/c [IS] 
and extends up to pt = 17 GeV/c. Preliminary results of this analysis have been 
published in [2H1 ED] • Also presented in this thesis arc the first measurements 
by STAR of 77 meson production. 
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1.4 Proton-proton collisions 

In QCD, the hadronic interactions are described in terms of the interactions of 
their constituent partons. The inclusive cross section of the reaction 

AB CX 

is calculated as the weighted sum of differential cross sections of all possible 
parton scatterings that can contribute [6]: 

^ d^a{AB-^CX) ^ f\ f\ „„, ... , 1 daiab^cd) x 

"^PC abcd-^O Jo ^^c at 

(1.12) 

Here f'^{xa) is the parton density function (PDF) that gives the probability 
that hadron A contains a parton a which carries the fraction Xa = qa/PA of 
its momentum. A similar definition applies to the density fg{xb)- The cross 
section da/dt of the hard partonic scattering 

ab cd 

is calculated in pQCD. The invariant kinematic variables for the partonic sub- 
process are 

s = {qa + qbf 
(qa- qcf, 

where is the partonic center-of-mass energy and is the momentum 

transfer from a to c. The fragmentation function D^'{z) in Eq. (|1.12p describes 
the probability that a given parton c produces a final state hadron C carrying 
a momentum fraction Zc = pc/qc- 

It follows from the above that the cross section calculations rely on two 
inputs — parton densities and fragmentation functions . These functions 
are non-perturbative, so that they cannot be calculated in QCD from first prin- 
ciples. However, they represent a properties of individual hadrons independent 
of the process in which they participate. Parton densities and fragmentation 
functions can, therefore, be obtained from an analysis of a large variety of scat- 
tering data. 

A widely used set of parton densities is obtained by the CTEQ Collabo- 
ration from a global QCD analysis of a large body of experimental data ;3lj . 
The global fit, together with a detailed treatment of published experimental 
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Figure 1.7: Nuclear shadowing effects on the R^a ratio, calculated with EKS98 |37) . nDS |38) . 
and HIJING [39] shadowing parametrizations. 



uncertainties, resulted in an excellent agreement with all available data. An 
alternative popular parametrization is MRST [32] . 

The fragmentation functions D'^{z) can be obtained directly from the pro- 
cess e+e" — > (7, Z) — > CX, in which the initial state has no hadrons. Such 
annihilation processes have been measured at many e+e~ colliders over a wide 
range of center-of-mass energies. The most recent parametrizations of fragmen- 
tation functions are KKP 33J, BKK [34 , BFGW 35], and Kretzer 36J. 

The cross sections for the individual partonic sub-processes are calculated 
in pQCD with no additional input, except for the strong coupling constant as. 
These calculations are usually performed at next-to-leading order (NLO), or 
even at next-to-next-to-leading order (NNLO). 

An important initial state effect in the heavy ion collisions is the modifica- 
tion of parton distribution functions inside nuclei. It is well known, that the 
quark structure functions at low fractional momentum are depleted in a nucleus 
relative to a free nucleon. This depletion is commonly referred to as nuclear 
shadowing. In Figure 11.71 we show the shadowing effects in d -I- Au collisions 
on the RdA ratio j40] , calculated with various parametrizations — EKS98 [37] , 
nDS [3B], and HIJING [32] ■ It is also a motivation for the present analysis to 
observe the nuclear shadowing and differentiate between models, although the 
required experimental precision may be prohibitively high. 
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Chapter 2 

The experiment 



2.1 RHIC accelerator complex 

The STAR experiment is located at the Brookhaven National Laboratory (BNL) 
on Long Island, USA. An important part of the physics program of the 
Laboratory is carried out at the Relativistic Heavy Ion Collider (RHIC). This 
is a multipurpose colliding beam facility [HI \TU\ , capable of accelerating protons, 
deuterons, and heavy ions over a broad energy range from the injection energy 
per nucleon of 10 GeV up to the top energy of 100 GeV for heavy ions and 
250 GeV for protons. 

The layout of the accelerator complex is shown in Figure I^TTl Heavy ions are 
accelerated in the Tandem Van de Graaff accelerator, the Booster, the Alternat- 
ing Gradient Synchrotron (AGS), and in the RHIC accelerator itself. The Linac 
serves to accelerate protons, which are then injected into the Booster. Below 
we will give a short description of each component of the accelerator complex. 

Tandem Van de Graaff generator Gold ions with unit negative charge 
are generated in the Pulsed Sputter Ion Source which delivers 250 jiA pulses 
of 600 IIS duration each. The ions are then accelerated in the Tandem Van de 
Graaff generator from the ground to +14 MV potential. They pass a set of 
stripping foils where they acquire a unit positive charge, and are subsequently 
accelerated again to the ground potential. The 1 MeV ions leaving the Tandem 
are stripped further to a charge of +32. There are two identical Tandems avail- 
able to provide two different ion species simultaneously (presently deuterium 
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Figure 2.1; Layout of the RHIC accelerator complex. Shown are the locations of the STAR, 
PHENIX, PHOBOS, and BRAHMS experiments around the RHIC ring. Also indicated in the 
figure is the STAR coordinate system with the positive z axis pointing in the West direction. 
Figure taken from [9] . 
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and copper in addition to gold). 

LINAC The LINAC serves to accelerate protons to an energy of 200 MeV, 
which are injected directly into the Booster. 

Booster synchrotron The 600 fis long Tandem pulse is injected into the 
Booster, after which the particles are captured into six bunches and accelerated 
to an energy of 95 MeV. Gold ions, when they are extracted from the Booster, 
are stripped to the charge +77, leaving only two tightly bound i^-shell electrons 
to be stripped at a later stage in the acceleration chain. 

AGS From the Booster, 24 bunches are injected into the AGS and rearranged 
into four final bunches containing 10^ ions each. Those bunches are accelerated 
to an energy of about 10 GcV. When transferred to the RHIC accelerator, the 
ions are fully stripped to the charge +39 in case of copper and to +79 in case 
of gold. 

RHIC accelerator The final stage of acceleration takes place in the RHIC 
synchrotron, where beams are circulating in two rings in opposite directions. 
The rings have a circumference of 3.83 km and are equipped with independent 
bending and focusing magnets and RF cavities. This provides the capability 
of operating the accelerator with two beams of unequal species. The bending 
magnets are superconductive and cooled by liquid helium. The complete cooling 
of the rings from room temperature to the operating temperature of 4.6°K takes 
about ten days. 

Up to 120 bunches can be injected in each ring and accelerated to an energy 
between 30 and 100 GeV. After acceleration, the bunches are transferred to the 
storage RF system, which maintains the bunch length at 1.52 m or 5 ns. The 
lifetime of a stored beam is about 10 hours, whereafter the beam is dumped 
and a new fill begins. A chosen pattern of empty buckets provides a sample of 
unpaired bunches crossing each interaction region for beam-background studies. 

Beams are made to cross at six points along the ring, four of which are used 
by the experiments STAR Jl], PHENIX [42^, PHOBOS [43^^, and BRAHMS [44]. 
Of the remaining two crossing points, one is occupied by the RF system, while 
the other is not used at present. 
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Table 2.1: RHIC runs in the years 2000-2007. 

Run Year Particle Beam energy Integrated Average beam 

species per nucleon luminosity polarization 











[GeV] 


[pb-i] 




[%] 


Run-1 


2000 


Au- 


h Au 


27.9 


< 10-9 










Au- 


h Au 


65.2 


20 X 10" 


-6 




Run-2 


2001-2002 


Au- 


h Au 


100.0 


258 X 10" 


-6 








Au- 


h Au 


9.8 


0.4 X 10" 


-6 








P- 


hp 


100.0 


1.4 




14 


Run-3 


2002-2003 


d- 


1- Au 


100.0 


73 X 10- 


-3 








P- 


\-p 


100.0 


5.5 




34 


Run-4 


2003-2004 


Au- 


h Au 


100.0 


3740 X 10" 


-6 








Au- 


h Au 


31.2 


67 X 10" 


-6 










I-P 


100.0 


7.1 




46 


Run- 5 


2004-2005 


Cu- 


hCu 


100.0 


42.1 X 10" 


-3 








Cu- 


hCu 


31.2 


1.5 X 10" 


-3 








Cu- 


hCu 


11.2 


0.02 X 10" 


-3 








P- 


hp 


100.0 


29.5 




46 






P- 


hp 


204.9 


0.1 




30 


Run-6 


2006 


P- 


hp 


100.0 


93.3 




58 






P- 


hp 


31.2 


1.05 




50 


Run-7 


2006-2007 


Au~ 


h Au 


100.0 


7250 X 10" 


-6 





RHIC performance To date, RHIC has delivered a variety of colliding beams 
of protons, deuterons, copper (Cu^^^), and gold ions (Au^^^) pT] [12]. 
In Table 12.11 we list the RHIC runs from the beginning of operations in the 
year 2000 up to the year 2007. The p + p runs provide reference data for the 
heavy ion physics program, as well as data to measure the proton spin structure 
at RHIC. For the latter purpose, the proton beams are polarized, reaching 
degrees of up to 60% in 2006. The data used in this thesis were taken in the 
d + Au run in 2002/03 and p + p in 2005, both at center-of-mass energies of 
200 GeV. 
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Silicon 




Figure 2.2: Perspective view of the STAR detector. Figure taken from I41| . 



2.2 STAR detector 

The STAR detector (Solenoidal Tracker At RHIC) [31] was designed primar- 
ily for measurements of hadron production in heavy ion and proton-proton 
colhsions over a large solid angle. For this purpose, large acceptance high 
granularity tracking detectors are placed inside a large volume magnetic field. 
A perspective view of the detector is shown in Figure 12.21 ^-nd a cutaway 
side view in Figure 12.31 

The barrel tracking detectors in STAR are a Silicon Vertex Tracker, sur- 
rounding the beam pipe, (SVT, not used in this analysis) and a large volume 
Time Projection Chamber (TPC), with an inner radius of 0.5 m, an outer ra- 
dius of 2 m, and a length of 4.2 m. The TPC covers a pseudorapidity range 
of I77I < 1.8 and is designed to reconstruct the very high multiplicity events 
produced in heavy ion collisions. These multiplicities can reach up to 1000 
charged tracks per unit rapidity in a central Au + Au collision at the largest 
beam energies. High granularity tracking in the forward and backward regions is 
achieved by two Forward TPCs (FTPC), each covering a range of 2.5 < [t;! < 4 
in pseudorapidity. 

For trigger purposes, the TPC is surrounded by a layer of scintillator tiles 
(Central Trigger Barrel, CTB, not used in this analysis) . 
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Magnet 




To trigger on the energy deposited by high transverse momentum photons, 
electrons, and electromagnetically decaying hadrons, a Barrel Electromagnetic 
Calorimeter (BEMC) [33] was incrementally added to the STAR setup from the 
year 2001 to 2005. 

The calorimeter surrounds and covers the full acceptance of the TPC and 
CTB. An Endcap Electromagnetic Calorimeter (EEMC) [M] was installed in 
2002-2003 to cover the pseudorapidity range 1 < 77 < 2. In the data taking 
period covered by this thesis, only the West half of the BEMC was fully opera- 
tional (0 < ?7 < 1). 

The STAR barrel detectors are placed inside a room temperature solenoidal 
magnet with maximum field of 0.5 T. The inner dimensions of the magnet are 
5.8 m in length and 5.27 m in diameter. 

To provide a minimum bias trigger and to measure centralities in heavy ion 
collisions, two sampling calorimeters (ZDC) are placed in the RHIC tunnel at 
18 m from the interaction point. Tiled arrays of scintillator counters (Beam- 
Beam Counter, BBC) are mounted around the beam pipe at a distance of 3.7 m 
from the interaction point, to provide a minimum bias trigger in p + p collisions. 
The detector subsystems relevant for the present analysis are briefly described 
in the following sections. We refer to Chapter [3] for a detailed description of the 
BEMC, which plays a central role in the analysis. 

Throughout this thesis we will use a Cartesian coordinate system defined 
as follows: z pointing along the beam in the West direction (see Figure [TT]) . 
y pointing upward, right-handed. 
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Figure 2.4: Schematic perspective view of the STAR TPC. Figure taken from 1471 . 

2.2.1 Time Projection Chamber 

The Time Projection Chamber [47 is the central tracking device in STAR. It 
allows one to track charged particles, measure their momenta, and identify the 
particle species by measuring the ionization energy loss dE/dx. 

A schematic layout of the TPC is shown in Figure 12.41 The TPC barrel 
measures 4.2 m in length and has an inner radius of 0.5 m and an outer radius of 
2 m. The TPC acceptance covers ±1.8 units in pseudorapidity and full azimuth. 
Particles are identified over a momentum range from 100 MeV/c to 1 GeV/c, and 
their momentum is measured in the range from 100 MeV/c to 30 GeV/c. 

The TPC is a gas filled cylindrical volume with a well defined uniform electric 
field gradient of about 135 V/cm. The secondary electrons released by ionizing 
particles along their path drift in the electric field towards the readout end- 
caps. The electric field is generated between a central membrane held at 28 kV 
potential and the endcaps, which are held at ground potential. A uniform field 
gradient is maintained by concentric equi-potential field cage cylinders biased 
via 2 Mfi resistors. The drift volume is filled with a gas mixture of 10% methane 
and 90% argon, which is held slightly above atmospheric pressure. The drift 
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FOIL WINDOW 

Figure 2.5: Perspective view of a STAR FTPC. Figure taken from |48| . 



velocity is 5.45 cm//is, and the maximum drift time ft'om the central membrane 
to endcap is 38.5 fis. 

The endcaps are instrumented with Multi-Wire Proportional Chambers 
(MWPC) with pad readout. The transverse coordinates of a track are recon- 
structed from the hits in the MWPCs, while the z coordinate is reconstructed 
from a measurement of the drift time. The total drift time of 5.45 fis is sampled 
by the readout electronics in 512 time buckets. 

In each endcap, the MWPCs are arranged in 12 sectors, each consisting of 
inner and outer sub-sector. The inner sub-sectors are in the region of highest 
track density and are, therefore, optimized for better two-track resolution, while 
the outer sub-sectors are optimized for better performance in the measurement 
of dE/dx. 

In the analysis presented in this thesis, the TPC is used as a charged par- 
ticle veto in the identification of photons in the BEMC. Samples of electrons 
reconstructed in the TPC serve to calibrate the energy response of the BEMC. 

2.2.2 Forward TPC modules 

Two Forward Time Projection Chambers (FTPC) [48] extend the STAR tracking 
capability to the pseudorapidity range 2.5 < jryj < 4. The layout of the 
FTPC is shown in Figure 12.51 Each FTPC is a cylindrical volume with a 
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Figure 2.6: Transverse view of the collision region, indicating the beam pipes, the Zero Degree 
Calorimeter, and the impact regions of neutrons and charged fragments with Z/A = 1 (top). 
Top view, showing the position of the ZDC modules behind the DX dipole magnets in between 
the two RHIC beam pipes (bottom). Figure taken from |49| . 



diameter of 75 cm and a length of 120 cm, with radial drift field and pad read- 
out chambers mounted on the outer cylindrical surface. Two such detectors are 
installed partially inside the main TPC on both sides of the interaction point. 
The FTPC is capable of reconstructing all charged tracks (typically 1000) travers- 
ing the detector in a central Au -|- Au event. 

In this thesis, the forward charged track multiplicity recorded in the FTPCs 
is used as a measure of the centrality in d -I- Au collisions. 

2.2.3 Zero Degree Calorimeter 

In addition to the STAR barrel detectors, a sampling calorimeter is placed at a 
distance of 18 m from the interaction point in the RHIC tunnel on both sides of 
the experimental hall, as shown in Figure BTBl These Zero Degree Calorimeters 
(ZDC) |501 149j are used to provide the minimum bias trigger and to measure 
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Figure 2.7: Mechanical design of the ZDC modules. Figure taken from I49| . 



centralities in heavy ion collisions. Furthermore, identical ZDC detectors are 
installed at each of the four RHIC experiments, providing comparable collision 
rate measurements to monitor the RHIC luminosity. 

The ZDC detector measures the total energy of the unbound neutrons emit- 
ted from the nuclear fragments after a collision. The charged fragments of the 
collision are bent away by the RHIC dipole magnets DX. In the upper plot of 
Figure [2?6l is shown a transverse view at the front face of the ZDC, indicating 
the position of the two beam pipes, the neutron spot inside the ZDC acceptance, 
and the spot of deflected fragments with Z/A = 1. 

The mechanical layout of the ZDC is shown in Figure 12.71 It consists of 
alternating layers of tungsten absorber and Cherenkov fibers with a total length 
of about 0.7 m. The transverse dimension oi x x y = 10 x 13.6 cm^ corresponds 
to an angular acceptance of about 2.5 mrad around the forward direction. 

In this thesis, we do not use the ZDC for centrality measurement, and refer 
to [IH] for details on such a measurement. For the d + Au data used in the 
present analysis, the ZDC provided a minimum bias trigger by requiring the 
detection of at least one neutron in the Au beam direction. The acceptance of 
this trigger corresponds to 95±3 % of the total d+An geometrical cross section, 
as determined from detailed simulations of the ZDC acceptance |TS]. 
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2.2.4 Beam-Beam Counter 

To provide a minimum bias trigger in p + p collisions, Beam-Beam Counters 
(BBC) [511 El] are mounted around the beam pipe beyond both poletips of the 
STAR magnet at a distance of 3.7 m from the interaction point. The BBC also 
serves to reject beam-gas events at the trigger level and to measure the beam 
luminosity in p -\- p runs. 

The detector consists of two sets of hexagonal scintillator tiles, see Figure l^TSl 
A ring with radius between 9.6 and 48 cm is formed by 18 small tiles, while 18 
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large tiles on the outside cover a radius between 38 and 193 cm. The small 
and large tile arrangements cover the pseudorapidities 3.4 < \r]\ < 5.0 and 
2.1 < \r]\ < 3.6, respectively. 

In p + p runs, a minimum bias trigger is provided by a coincidence of signals 
in at least one of the 18 small BBC tiles on each side of the interaction region. 

The two BBC counters also record the time of flight, which provides a mea- 
surement of the z position of the interaction vertex to an accuracy of about 
40 cm. Large values of the time of flight difference between the two BBC coun- 
ters indicate the passage of beam halo, which is rejected at the trigger level. 

A measurement of the counting rate in the BBCs allows for a determination 
of the absolute luminosity to an accuracy of about 15%, the relative luminosities 
per run are determined to a precision of better than lO'^ [HI [5^ . 



Chapter 3 



STAR Electromagnetic 
Calorimeter 

The Barrel Electromagnetic Calorimeter (BEMC) ^45j is a lead-scintillator 
sampling calorimeter, surrounding the STAR TPC as shown in Figure [3Tl The 
BEMC was installed in several stages during the period of 2001-2005. Only the 
West half of the BEMC was fully operational during the 2003 and 2005 runs 
which provided the data presented in this thesis. The Endcap Calorimeter |46j . 
which is not used in the present analysis, was installed in the years 2002-2003. 

The BEMC is used to trigger on and to measure jets, leading hadrons, 
direct photons, and electrons from heavy quarks produced at large transverse 
momentum. For this purpose, the BEMC provides large acceptance for photons, 
electrons, tt^, and rj mesons in all colliding systems ranging from p + p up to 
Au + Au. In the next sections we will describe the BEMC in more detail. 

3.1 Mechanical layout 

The calorimeter is located inside the magnet coil and surrounds the TPC. It 
covers a pseudorapidity range of |ry| < 1 and full azimuth, matching the TPC 
acceptance. The calorimeter is divided in two adjacent barrels, one positioned 
at the West half of the STAR detector (0 < 77 < 1) and the other one at the 
East half (— 1 < 77 < 0). Each half-barrel has a length of 293 cm, an inner radius 
of 223 cm, and an outer radius of 263 cm. 

The half-barrel is azimuthally segmented into 60 modules. Each module is 
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Figure 3.1: Cross-sectional and longitudinal view of the STAR detector, showing the layout 
of the BEMC. Figure taken from |45) . 



approximately 26 cm wide and covers 6 degrees (17 mrad) in azimuth and one 
unit in pseudorapidity. The active depth is 23.5 cm, to which is added 6.6 cm 
of structural elements at the outer radius. The longitudinal and transverse 
segmentation of a module is shown in Figure 13.21 and the radial structure in 
Figure 13.31 

The modules are segmented into 40 projective towers of lead-scintillator 
stacks, 2 in (/) and 20 in rj. A tower covers 0.05 in A(p and 0.05 in Arj. Each 
calorimeter half is thus segmented into a total of 2400 towers. 

Each tower consists of an inner stack of 5 layers of lead and 5 layers of 
scintillator, and an outer stack of 15 layers of lead and 16 layers of scintillator. 
All these layers are 5 mm thick, except the innermost two scintillator layers, 
which are 6 mm thick. A separate readout of these latter two layers provides the 
calorimeter preshower signal. A Shower Maximum Detector (SMD) is positioned 
between the inner and outer stacks, at a depth of appoximately 5 radiation 
lengths. The whole stack is held together by mechanical compression and friction 
between layers. 




Figure 3.2: Side view of a calorimeter module and top view of a scintillator plate segmented 
into 20 X 2 towers. Figure taken from |45) . 
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Figure 3.3: Transverse view of a calorimeter module, showing the inner layer of 
lead/scintillator stacks, the shower maximum detector (SMD), the outer layer of stacks, and 
the carriage structure at the outer radius of the BEMC barrel. Figure taken from 1451 . 
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Figure 3.4: Optical readout scheme of a BEMC tower. Figure taken from I45| . 



3.2 Optical structure 

The plastic scintillator layers are machined as "megatiles" , covering the full 
length and width of a module. These megatiles are segmented into 40 optically 
isolated tiles, as shown in the top diagram of Figure [3?2l The optical separation 
between the individual tiles is achieved by 95% deep cuts in the scintillator filled 
with opaque epoxy. The optical crosstalk between adjacent tiles is reduced to 
a level of 0.5% by painting a black line on the surface opposite to the isolation 
groove. 

The optical readout scheme is shown in Figure 13.41 The signal from each 
tile is collected by a wavelength shifting (WLS) fiber embedded in a cr-groove in 
the tile. The WLS fibers run along the outer surface of the stack and terminate 
in an optical connector mounted at the back-plate of the module. From the 
back-plate, 2.1 ni long fibers run through the STAR magnet structure to the 
readout boxes mounted on the outer side of the magnet. In these boxes, the 21 
fibers from the tiles of one tower are connected to a single photomultiplier tube 
(PMT). The PMTs are powered by Cockroft-Walton bases, which are remotely 
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controlled over a serial communication line by the slow control software. 

From layer by layer tests of the BEMC optical system, together with an 
analysis of cosmic ray and test beam data, the nominal energy resolution of the 
calorimeter is estimated to be SE/E = 15%/^£:[GeV] © 1.5% [S3]. 

3.3 Shower Maximum Detector 

The Shower Maximum Detector (SMD) is a multi-wire proportional counter 
with strip readout. It is located at a depth of approximately 5.6 radiation 
lengths at ?7 = increasing to 7.9 radiation lengths at 77 = 1, including all 
material immediately in front of the calorimeter. 

The purpose of the SMD is to improve the spatial resolution of the calorime- 
ter. This is necessary because the transverse dimension of each tower (about 
10 X 10 cm^) is much larger than the lateral spread of an electromagnetic shower. 
The improved resolution is essential to separate the two photon showers origi- 
nating from the decay of high momentum tt'^ and ry mesons. 

The layout of the SMD is shown in Figure 13.51 Independent cathode planes 
with strips along 77 and directions allow the reconstruction of a two-dimensional 
image of a shower. The coverage in Arj x is 0.0064 x 0.1 for the 77 strips 
and 0.1 x 0.0064 for the strips. There are a total of 36000 strips in the full 
detector. 

Beam test results at the AGS have shown that the SMD has an approxi- 
mately linear response versus energy. The energy resolution in the 77 coordinate 
(front plane) is approximately 5E/E = 86%/-\/E[GeVj 12%, whereas that in 
the 4> coordinate (back plane) is worse by about 3-4%. The position resolution 
is a{r(j)) = 5.6/a/£;[GcV] e2.4mm and cr(z) = 5.8/A/E[GeVj © 3.2 mm. 

3.4 Preshower Detector 

The first and second scintillating layers of each calorimeter module are used 
as a preshower detector (PSD). To achieve a separate readout of these layers, 
two WLS fibers are embedded instead of one in the cr-groove of each tile. This 
additional pair of fibers from the two layers illuminate a single pixel of a multi- 
anode PMT. A total of 300 16-pixel multi-anode PMTs are used to provide the 
4800 tower preshower signals. 

The preshower detector was fully instrumented and read out only in 2006, 
so that it could not be used in the present analysis. 
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Figure 3.5: Schematic illustration of the SMD, showing in the top figure a three dimensional 
view of the extruded aluminium profile, containing the anode wires, and the two readout 
pad planes, running parralel and perpendicular to the wires. The profile of a BEMC shower, 
as recorded in the two SMD pad planes, is shown by the histograms. The bottom plot 
shows a transverse view of the aluminium extrusion, the anode wires, and the pad planes. 
Figure taken from |45| . 
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3.5 BEMC electronics 

The calorimeter is a "fast" detector in STAR, so that its ADCs can be read 
out on each RHIC bunch crossing. The calorimeter data is also used in the 
STAR level-0 trigger, in the form of the "High Tower" and "Patch Sum" trigger 
primitives. 

The level-0 HighTower trigger used in this analysis is a requirement that 
the energy deposited in any single calorimeter cell in the event exceeds a given 
threshold. This allows one to enhance the statistics at the high energy part of 
the spectrum. 

The complete description of the BEMC electronics operation is given in 
Appendix |^ 



Chapter 4 

Event reconstruction in 
STAR 

4.1 Data aquisition and trigger 

The STAR data aquisition system (DAQ) [54] receives the input from multiple 
detectors at various readout rates. The typical recorded event rate of 100 Hz is 
limited by the drift time in the TPC (the slowest detector in STAR). The total 
event size can reach up to 200 MB in Au + Au collisions. STAR takes data in 
runs of about half an hour duration, each having 50-100 x 10"^ events. 

The STAR trigger [5S] is a pipelined system, capable to cope with the RHIC 
beam crossing frequency of 10 MHz. The trigger processes information from 
fast detectors, such as the ZDC, BBC, CTB, or BEMC, and decides if the event 
should be read out and saved to tape. Each event is categorized by multiple 
trigger criteria, and the events selected by various branches of the decision tree 
are written to tape, sharing the available 100 Hz DAQ bandwidth. 

The datasets used in the present analysis were taken in the d+Au run of 2003 
and the p + p run of 2005, see also Table [2TT1 The following trigger conditions 
had to be satisfied: 

Minimum bias (MinBias) trigger in d + Au collisions This condition 
required the presence of at least one neutron signal in the ZDC in the gold beam 
direction. As given in [T^ , this trigger condition captured 95 ± 3 % of the total 
d + Au geometric cross section of CThack^" = 2.21 ± 0.09 b. 
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Table 4.1: HighTower trigger thresholds used in d + Au 2003 and p + p 2005 
data runs. 

Dataset HighTower-1 threshold HighTower-2 threshold 

[GeV] [GeV] 

d + Au 2003 2.5 4.5 

p + p 2005 2.6 3.5 



MinBias trigger in p+p collisions This condition required the coincidence 
of signals from two BBC tiles on the opposite sides of the interaction point. 
Due to the dual-arm configuration, this trigger is sensitive to the non-singly 
diffractive (NSD) cross section, which is a sum of the non-diffractive and doubly 
diffractive cross section. The total inelastic cross section is a sum of the NSD 
and singly diffractive cross section. 

A minimum bias cross section of ctbbc = 26.1 ± 0.2(stat) ± 1.8(syst) mb was 
independently measured via Vernier scans in dedicated accelerator runs [5S] . 
This trigger captured 87 ± 8 % of the p + p non-singly diffractive (NSD) cross 
section, as was determined from the detailed simulation of the BBC accep- 
tance [T^. Correcting the BBC cross section for the acceptance, we obtain a 
value for the NSD cross section of ct^sd = ^^^-^ ^ '^■^ 

HighTower trigger This condition required, in addition to the MinBias, an 
energy deposit above a predefined threshold in at least one calorimeter tower. 
The purpose of this trigger is to enrich the sample with events that have a large 
transverse energy deposit. Two different thresholds were applied, giving the 
HighTower-1 and HighTower-2 datasets. The values of these thresholds for the 
various runs are shown in Table [4A] 



4.2 STAR reconstruction chain 

The events recorded on tape are passed through the standard STAR reconstruc- 
tion chain. This reconstruction is performed routinely on the RHIC Computing 
Facility (RCF), which is a large computing farm located at BNL. 

The most important part of the data reconstruction at this stage is track- 
ing in the TPC and FTPCs. Charged tracks are reconstructed in the main 
TPC using a Kalman filter [57j . and in the FTPCs using a conformal mapping 
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method [5S] . The primary vertex is found by extrapolating and intersecting all 
reconstructed tracks. The vertex resolution in z is between 0.9 mm and 0.35 mm 
depending on the track multiplicity, whereas in the transverse plane it is about 
0.5 mm. Once the vertex has been found, all tracks that approach to it closer 
than 3 cm are re-fitted to include the vertex position as the origin. Although the 
wire chambers are sensitive to almost 100% of the drifting electrons, the overall 
tracking efficiency is only 80-90% due to fiducial cuts, track merging, bad pads, 
and dead channels. The momentum resolution of tracks worsens linearly with 
PT from 2% for 300 MeV/c pious to 7% for 4 GeV/c pious. 

Because the BEMC reconstruction is not yet performed in the standard 
STAR reconstruction, the raw BEMC data are passed to the physics analysis. 
The tower ADC data can be directly passed because they only take a small 
fraction of the event size. The SMD strip ADCs are zero-suppressed and then 
also passed to the analysis. This scheme implies that removal of malfunctioning 
elements and a full calibration of the BEMC is performed as a part of the physics 
analysis. This has the advantages that the reconstructed electron tracks in the 
TPC can be used to calibrate the energy response of the BEMC, and that the 
successive improvements in the BEMC calibration do not require re-generating 
the full dataset from the raw events on tape. 

All data reconstruction and analysis in STAR is performed using the ROOT 
framework [55]. The processing of a full dataset, such as p + p oi d + Au, takes 
about three months. 

4.3 BEMC status tables 

A quality assurance (QA) procedure for the BEMC is routinely performed 
before the physics analysis, in order to remove malfunctioning detector com- 
ponents from the data and to correctly reproduce the time dependence of the 
detector acceptance in the Monte Carlo simulation. This QA procedure results 
in timestamped status tables, which are used as an input to the physics analysis. 
Below we describe the QA procedure performed for the BEMC towers, a similar 
procedure is apphed to the SMD strips. 

For each run, the raw ADC spectra of all towers were accumulated and a 
number of criteria were applied to recognize common failure modes, such as the 
malfunctioning of entire readout boards and crates. A typical ADC spectrum 
of a tower is displayed in Figure 14.11 and shows the signal distribution and the 
accumulation of ADC counts in absense of a signal (pedestal). The position 
of these pedestals provide the zero offset of the ADC measurement and are. 



38 



CHAPTER 4. EVENT RECONSTRUCTION IN STAR 




Figure 4.1: ADC spectrum of a BEMC tower, showing a pedestal located at 30 ADC counts. 



together with the width, stored in time dependent tables for each tower. Chan- 
nels with anomalous pedestal positions and widths are flagged as bad. The 
signal fraction was defined as the number of ADC counts that are more than six 
standard deviations above the pedestal. Towers with a signal fraction smaller 
than 0.001 are flagged as "cold" or "dead" , while those with a fraction above 0.1 
are marked as "hot" or "noisy" (the exact numbers are multiplicity dependent 
and are adjusted for each collision system). Saved are, as function of run num- 
ber, the position of the pedestals, their widths, and flags indicating the status 
of each tower. The average fraction of good towers was found to be about 90% 
in the 2003 d4- Au run, with run-to-run fluctuations of about 2-5%. In the 2005 
p + p data the fraction of good towers was found to be about 97%. 

In the BEMC reconstruction performed in this analysis, the status tables 
were read in and used for pedestal subtraction of the ADC signals and for 
removal of towers which were flagged as bad. 

4.4 BEMC energy calibration 

The purpose of the energy calibration is to establish the relation between ADC 
counts and the energy scale in GeV. The calibration proceeds in two stages. 
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First, a relative calibration matches the gains of individual towers to achieve 
an overall uniform response of the detector. A common scale between ADC 
counts and energy is then determined in a second absolute calibration step. 
The relative tower-by-tower calibration is done using minimum ionizing particles 
(MIP), while the absolute energy scale is determined from energy measurements 
of identified electrons in the TPC. 

4.4.1 MIP calibration 

A significant fraction (30-40%) of high energy charged hadrons traversing the 
BEMC only deposit a small amount of energy in the towers, equivalent to a 
250-350 MeV electron, largely due to ionization energy loss (minimum ionizing 
particles). The signal from these particles is usually well separated from the 
tower pedestals. 

To identify MIP particles, TPC tracks of sufficiently large momentum above 
1.2 GeV/c are extrapolated to the BEMC and the response spectra are accumu- 
lated, provided that the track extrapolation is contained within one tower and 
that there are no other tracks found in a 3 x 3 patch around this tower. In 
Figure [321 is shown a tower ADC spectrum collected from the d + Au dataset, 
which clearly shows the position of the MIP peak superimposed on a broad 




Figure 4.2: The BEMC tower response to MIP tracks. The curves indicate a fit to two 
Gaussians, one for the peak and one for the background. Figure taken from [60) . 
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Ptpc |GeV/c) 

Figure 4.3: Electron energy measured in BEMC vs. momentum measured in TPC. The first 
order polynomial fit determines the global calibration constant. Figure taken from |60j . 

background [50] . The position of the fitted Gaussian is calculated for each tower 
and used to calculate the tower-by-tower gain corrections needed to equalize the 
detector. 

The disadvantage of this method is that the calibration is performed at the 
low end of the scale, where the signal is more susceptible to noise and where the 
lack of lever arm does not allow to detect possible non-linearities in the detector 
response. 

4.4.2 Electron calibration 

Because the electron momentum can be independently measured in the TPC, 
it is possible to calibrate the absolute energy scale of the calorimeter using the 
simple relation for the ultra-relativistic electrons, E/p = 1. 

Figure 14.31 shows the electron energy measured in the calorimeter versus 
its momentum measured in the TPC [60 . The calorimeter response is quite 
linear up to 8 GeV, and the global gain correction obtained from the linear fit 
is applied to all towers. 

This method takes advantage of the well understood TPC detector for the 
precise measurement of the electron track momentum in a wide range. However, 
it requires high statistics to calibrate the high energy part of the spectrum, 
so that only one global calibration constant for the calorimeter is obtained at 
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present. The systematic year-to-year uncertainty on the electron calibration 
was estimated to be 5% [5T] . 

It has been found that the current calibration is less reliable at the edges of 
the calorimeter half-barrel, therefore, the tower signals from the two 77-rings at 
each side are later removed from this analysis (see Section !??^ . 

This combination of the MlP-based equalization and electron-based absolute 
calibration is applied to the data after each running period, starting from 2003 
d + Au run. The run dependent calibration constants are saved in the STAR 
database and automatically applied to the ADC readout in the software. 

4.5 Event selection 

The event selection starts with rejecting events where subdetectors needed for 
this analysis were not operational or malfunctioning. In the following sub- 
sections we will describe several additional selection criteria in detail. 

4.5.1 Beam background rejection 

In d + Au events, interactions of gold beam particles with material approxi- 
mately 40 m upstream from the interaction region give rise to charged tracks 
that traverse the detector almost parallel to the beam direction. To identify 
events containing such background tracks, the ratio 

-Ebemc 
r = 

-Bbemc + Etpc 

is calculated, where Ebkmc is the total energy recorded in the BEMC and i?TPC 
is the energy of all charged tracks reconstructed in TPC. In events containing 
background tracks, the ratio r tends to become large because the background 
tracks give a large energy deposit in a calorimeter without being reconstructed 
in the TPC, since they do not point to the vertex. This is shown in Figure [131 
where the distribution of r is plotted for the d-l- Au and p+p datasets. The peak 
near unity in the left-hand plot indicates the presence of beam halo in d + Au 
collisions, and events with r > 0.8 were removed from the d + Au analysis. This 
cut rejected 3.4% of MinBias and 13% of HighTower-2 triggered events. From 
a polynomial fit to the d + Au distribution in the region r = 0.6-0.8 (curve in 
Figure the false rejection rate was estimated to be 3.6% in the d + Au 

HighTower-2 data and less than 1% in the other datasets. 
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■^BEMC / (■'^BEMC ' -^TPc) -^BEMC / (■'^BEMC ' -^TPc) 



Figure 4.4: Distribution of r = i?BEMc/(^BEMC + £'TPc): which shows beam background at 
r > 0.8 in d+Au events (left) and its absense in p+p events (right). The curve in the left hand 
plot indicates a polynomial fit used to estimate the false rejection rate in the d + Au data. 



The cut was not applied to the p + p data since here the beam background 
is almost absent, as can be seen in the right-hand plot of Figure 

During the summer in 2006, additional shielding walls were installed in STAR 
to reduce this beam background to a negligible level. 

4.5.2 Vertex reconstruction 

The event vertex is reconstructed to an accuracy of better than a millimiter in 
the z direction, from the tracks reconstructed in the TPC. The distribution of 
the vertex z coordinate in the p + p MinBias data is shown in Figure 14.51 

Events with l^vortexl > 60 cm were rejected in the analysis, as indicated 
by the vertical lines in Figure 14.51 This cut is applied because the amount of 
material traversed by a particle increases dramatically at large values of l^vcrtcxl- 
As a consequence, the TPC tracking efficiency drops for vertices located far from 
the center of the detector. 

In the HighTower trigger data, the track multiplicity is almost always suf- 
ficient for a TPC vertex reconstruction, but this is not so in the p + p and 
d + Au minimum bias data. Since the p + p minimum bias trigger is based 
on coincidences in the BBC, we can use the timing information of the BBC to 
reconstruct a vertex for every event, even when the TPC vertex reconstruction 
failed (about 35% of the minimum bias events). The timing information from 
the BBC was calibrated against the z vertex coordinate reconstructed in the 
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Figure 4.5: Vertex distribution in the p + p MinBias dataset. Vertical lines indicate the cut 
used in the analysis. 



TPC, as illustrated in Figure [4!6l (top) . where we show the correlation between 
the BBC time difference At and Zvcrtox in the TPC. The straight line in the 
plot corresponds to a linear fit 

^vertex — aAi + b 

yielding a = 2.824 ± 0.003 cm per ADC count and b ^ 11.00 ± 0.02 cm. In the 
bottom plot of Figure [4?6l we show the distribution of zbbc ~ -^tpCi together 
with a Gaussian fit. From this fit we obtain the BBC vertex resolution of 40 cm. 

Whereas p + p events without a TPC vertex can be recovered by using the 
BBC timing information, this cannot be done for d + Au events because the 
BBC is not in the trigger and timing information may be absent. Since the 
7r° reconstruction requires the presence of vertex, the d + Au events without a 
TPC vertex are removed from the analysis. The vertex finding efficiency was 
determined from detailed Monte Carlo simulation of the full d + Au events and 
was found to be 93±1 % in the ±60 cm window [Kj. This result is used to correct 
the d + Au data for vertex inefficiencies, as will be explained in Section [7. II 
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4.5.3 HighTower trigger condition 

The HighTower-triggered data are filtered using a software impleraentation of 
the HighTower trigger. In this filter, the highest tower ADC value found in the 
event is required to exceed the same HighTower-1 (HighTower-2) threshold as 
the one that was used during the run. This filter is needed to remove events that 
were falsely triggered due to the presence of noisy channels (hot towers) . Such 
channels are identified offline in a separate analysis and recorded in a database 
as described Section 14.31 This software filter also serves to make the trigger 
efficiency for Monte Carlo and real data as close as possible. 

4.6 Centrality selection in d + Au data 

To measure the centrality in d + Au collisions, we use the correlation between 
the impact parameter of the collision and the charged track multiplicity in the 
forward direction. This correlation was established from a Monte Carlo Glauber 
simulation [Tni El] using, as an input, the Woods-Saxon nuclear matter den- 
sity for the gold ion [53] and the Hulthcn wave function of the deuteron [Mj . 
In this simulation, the inelastic cross section of an individual nucleon-nucleon 
collision was taken to be (t^^ ~ 42 mb. The produced particles were then 
propagated through a full GEANT simulation of the STAR detector and the 
charged track multiplicity was recorded, together with the number of nucleon- 
nucleon collisions simulated by the event generator. 

For the event-by-event centrality determination, we measured the multiplic- 
ity (A^FTPc) of tracks reconstructed in the FTPC-East acceptance (in the Au 
beam direction), following the centrality binning scheme used in other STAR 
publications [13 [SS]. The following quality cuts were applied to the recon- 
structed tracks: (i) at least 6 hits are required on the track; (ii) pT < 3 GeV/c, 
to guarantee that the track is fully contained in the FTPC acceptance, and 
(iii) distance of closest approach (DCA) to the vertex should be less than 3 cm. 
The multiplicity distributions obtained from the d + Au dataset are shown in 
Figure [4771 for the MinBias, HighTower-l , and HighTower-2 triggers. 

Based on the measured multiplicity, the events were separated into three 
centrality classes: 0-20% most central, 20-40% mid central, and 40-100% most 
peripheral, as illustrated by the vertical lines in Figure H771 

Table [121 lists the iVpTPC ranges that defined the centrality classes, and the 
corresponding mean number of binary collisions (A^cou) in each class, obtained 
from the Glauber model. In the table are also listed the systematic uncertain- 
ties on (iVcoii), which are estimated by varying the Glauber model parameters. 
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Figure 4.7: Centrality selection in d + Au data, based on the FTPC multiplicity. Three 
centrality classes are defined, containing 0-20% most central, 20-40% mid central, and 
40-100% most peripheral events. 



Table 4.2: Centrality classes defined for the d + Au data and the corresponding 
(.^coii) values [551 mi- The errors given for (iVcou) indicate the systematic 
uncertainty. 



Centrality class -/VpTPc range (-/Vcoii) 



d + Au minimum bias 






7.5 ±0.4 


0-20% most central 


< 


10 


15.0±1.1 


20-40% mid central 


10 


-16 


10.2±1.0 


40-100% most peripheral 


> 


17 


4.0±0.3 



p + p — 1 



Chapter 5 



Neutral meson 
reconstruction 

The goal of this analysis is to measure tt*^ and ij production in d + Au and p + p 
collisions. The 7r° and 77 are identified by their decay 

TT*^ — > 'fj and 1] — s- 77. 

These decay modes have branching ratios of 0.988 and 0.392, respectively [66] , 
The BEMC is used to detect the decay photons, as will be described in the next 
sections. The lifetime of the tt" is t = 8.4 x 10~^^ s, which corresponds to a 
decay length cr = 0.025 p,m. The lifetime of the 77 is even shorter (7 x 10~^^ s). 
Therefore, we can assume that the decay photons originate from the primary 
vertex. For each event, the invariant mass 

M^^ = ^2EiE2{l-cos'tp) (5.1) 

is calculated for all pairs of photons detected in the BEMC. Here Ei and E2 
are the energies of the decay photons and ^ is the opening angle between them, 
as measured in the laboratory system. 

The reconstructed masses are accumulated in invariant mass spectra, where 
the TT*^ and the ij show up as peaks around their nominal masses. These peaks 
are superimposed on a broad distribution of combinatorial background, which 
originates from photon pairs that are not produced by the decay of a single 
parent particle. 

In Table 15.11 we list the number of events in all datasets used in the analysis 
after the event selection procedures described in section HTSl were applied. 
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Table 5.1; Statistics used in the analysis after the event selection. 

Dataset Number of events 

MinBias HighTower-1 HighTower-2 
d + Au 164 608 53 154 40 974 

0-20% most central 21 382 12 567 8 744 

40-100% most peripheral 108 904 33 201 24 658 

pTp 4 433 817 920 567 872 811 



5.1 BEMC clustering 

The first step in the invariant mass reconstruction is to find clusters of energy 
deposits in the calorimeter. The purpose of the cluster finding algorithm is 
to group adjacent hits that are likely to have originated from a single incident 
photon. The algorithm is applied to the BEMC tower and preshower signals, 
as well as to the signals from each of the two SMD layers. 

The clustering algorithm starts by accumulating a list of cluster seeds that 
contains all hits in a module with an energy deposit above a certain threshold 
(i?sccd)- Starting from the most energetic seed in the list, an energy ordered list 
of module hits is searched for those adjacent to the present cluster. When such 
a hit is found, then, provided that it is above a second threshold (i?add), it is 
added to the cluster and removed from the list. The clustering stops when either 
a pre-defined maximum cluster size (iVmax) is reached or no more adjacent hits 
are found. The clustering algorithm then proceeds to process the next most 
energetic seed. At the end, clusters with total energy below the third threshold 
(i'min) are discarded. Note, that, by construction, the clusters are confined 
within a module and cannot be shared by adjacent modules. However, the 
likelihood of cluster sharing between modules is considered to be low since the 
modules are physically separated by about 12 mm air gaps. In Table we list 
the threshold values used in the clustering algorithm for all four detectors. 

In Figure 15.11 we show the assignments made by the algorithm on several 
possible one-dimensional cluster topologies. Note, that the rightmost hit pattern 
in this figure shows a double-peak structure, which is splitted into two adjacent 
clusters by the algorithm. However, statistical fluctuations in single photon 
signals may also be the cause of a double-peak structure. In such a case, the 
cluster splitting by the algorithm becomes a source of background, as will be 
discussed in Section [531 
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Table 5.2: Cluster finder threshold values used in the analysis. 



Detector 


-^sccd 


-Eadd 




max 




[GeV] 


[GeV] 


[GeV] 




Towers 


0.35 


0.035 


0.02 


4 


Preshower 


0.35 


0.035 


0.02 


4 


SMD-77 


0.2 


0.0005 


0.1 


5 


SMD-</) 


0.2 


0.0005 


0.1 


5 



I Not part of a cluster 
mi Cluster hits 

Order of association 
(sorted by energy) 




Figure 5.1: Schematic view of hit topologies in the BEMC and the assignment of hits to 
clusters by the algorithm described in the text. 
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The readout of the SMD 77 and 4> planes is one-dimensional, so that there 
is no ambiguity in what is considered to be an adjacent hit. The calorimeter 
tower readout is two-dimensional, and two hits are considered to be adjacent 
when they share a side and not when they share only a corner. 

The cluster position in the 77 and (p coordinates is calculated as the energy 
weighted mean position of the participating hits. In this calculation, the geo- 
metrical center of the detector element is taken as the hit position. 

After the tower, preshower, and SMD clusters are found, the next step is 
to combine them into so-called BEMC points that should correspond as closely 
as possible to the impact point and energy deposit of a photon that traversed 
the calorimeter. This procedure treats 2x2 tower patches corresponding to the 
SMD-(/) segmentation, as shown by the top diagram in Figure 15.21 

It is required that every reconstructed BEMC point contains a tower cluster, 
since the energy deposit of the incident particle is measured in the BEMC towers. 
Adding information from the SMD leads to a variety of combinations, as shown 
schematically by the diagrams (a)-(c) in Figure [5?2l In the following paragraphs 
we describe how each case leads to the reconstruction of a BEMC point. 



Tower, SMD-r; and SMD-</) clusters The algorithm calculates for all 
combinations of SMD-77 and SMD-0 clusters in a patch the energy asymmetry 

A^\E^~ E^\/{E^ + E^), 

where E^i and E^ are, respectively, the energy deposits measured in the SMD-7y 
and SMD-(/) planes. 

The cluster assignment constitutes a well known problem in combinatorics 
(Assignment problem [ST]) which we solve by a call to the CERN library routine 
ASSNDX [68] that combines objects into pairs in a way that minimizes the total 
cost. In the present algorithm the cost function is defined as energy asymmetry 
A between clusters. 

Each associated SMD pair is matched to the tower cluster closest in rj and 
(j). The total tower energy in a patch (including unassociated) is shared between 
points, weighted by their average SMD energy, that is, each i-th pair will produce 
a point with energy 

passoc 



l^j -^SMD, j 



where E^^ - = (i?,,, i + E^_ i)/2. The 77 and (j) coordinates are that of the 
SMD clusters. 
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10 patches per module, 2x2 towers 
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Figure 5.2: Combining tower and SMD clusters into BEMC points. 
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This procedure works well, provided that the occupancies of the 2x2 tower 
patches are low. Indeed, the number of tower or SMD clusters reconstructed 
even in the most central d + Au events is below 30 in the complete half- 
barrel, corresponding to a mean number of 12 clusters per event and an average 
occupancy of 2% per patch. 

Tower and SMD-?7 clusters In this case, the tower and SMD-ry clusters 
are associated by the same algorithm as used above, except that here the cost 
function is defined by the energy asymmetry 

A = |^t-^^|/(^t + -B^), 

where Et is the energy deposit in a tower and E^j is the energy deposit measured 

in the SMD-77 plane. The total energy of tower clusters in a patch is shared 
between associated pairs, weighted by their tower energy: 

^assoc 

Ei = £;J°**' X . 

t passoc 

The r] coordinate associated to the BEMC point is taken directly from the 
SMD-ry cluster, while the </> coordinate is taken from the tower cluster. 

Tower and SMD-<^ clusters This case is treated as described above. The 
resulting BEMC points will have the (j) coordinate from the SMD-t^ clusters and 
the T} coordinate from the tower clusters. 

Tower clusters only If there are no SMD clusters in a patch that contains 
the tower cluster position, the energy and coordinates of the BEMC point are 
taken to be those of the tower cluster. 

The relative occuranccs of these four cases arc approximately in proportion of 
70 : 5 : 5 : 20 % for clusters with energy above 4 GeV, and 25 : 10 : 10 : 55 % at 
the lower energies. 

The information about the shower shape in the SMD is in principle available 
but not used in the present clustering algorithm. 

5.2 BEMC cluster cuts 

After clustering, only the BEMC points containing tower and both SMD-t/ and 
SMD-(/) clusters were kept to be used in the further analysis of the HighTower- 
triggered data. In the analysis of MinBias data all reconstructed BEMC points 
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were used, even when they do not contain SMD clusters. From the decay kine- 
matics in the laboratory it follows that the opening angle between the photons is 
smallest when these photons equally share the energy of the parent. In Figure [5T51 
is shown this minimal opening angle versus the energy of the parent 7r° or i] and 
compared to the tower and SMD strip size. It is seen that the spatial resolution 
of better than a calorimeter tower is needed to resolve the decay photons of 
neutral pions with momenta larger than 5 GeV/c. For this reason, the SMD 
information is essential in this analysis. 

It is seen from beam tests |53j that the SMD efficiency decreases rapidly 
with energy of the traversing particle and is smaller than 50% at E < 2 GeV/c. 
The energy resolution SE/E = 12% © 86%/Ve is also poor at low energy, so 
that significant fluctuations in the strip readout are expected. Therefore, an 
SMD cluster is required to contain signals from at least two strips in order to 
be accepted in the HighTower-l data. This cut rejects a large fraction of the 
distorted and falsely split SMD clusters, and reduces a possible effect of poor 
SMD response simulation at low energies. 

It has been found that the tower calibration is less reliable at the edges of 
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the calorimeter acceptance. For this reason, we only keep the reconstructed 
clusters in the range 0.1 < rj < 0.9 for the further analysis, excluding two tower 
77-rings at each side of the calorimeter half-barrel. 

A charged particle veto (CPV) cut is applied to reject the charged hadrons 
that are detected in the calorimeter. These charged hadrons can be recognized 
as BEMC clusters with a pointing TPC track. The cluster was rejected if the 
distance between the BEMC point and the closest TPC track (D) was smaller 
than 0.04 in the r]-(p coordinates. 



D = v/A?72 + A02 < 0.04. 

The BEMC points remaining after this cut are considered to be photon candi- 
dates, which are combined into pairs, defining the set of vr" candidates. 

The asymmetry of the two-body decay of neutral mesons is defined as 

El + E2 

where Ei and E2 are the energies of the decay photons. From the decay 
kinematics it follows that this energy asymmetry is uniformly distributed 
between and 1 In Figure we show the distribution of the asymmetry of 
photon pairs reconstructed in p + p data. In the MinBias data the distribution 
is not flat because of the acceptance effects — photons from the asymmetric 
decay have a large opening angle and there is a large probability that one of 
them escapes the barrel. It is also seen that the HighTower energy threshold 
biases the asymmetry to the higher values, because it is easier for an asym- 
metric decay to pass the trigger. In this analysis, the 7r° candidates were only 
accepted if the asymmetry was less than 0.7, in order to reject very asymmetric 
decays, where one of the BEMC points has low energy, and to reject a signifi- 
cant part of the low mass background (this background will be described in the 
following sections) . It turns out that the asymmetry cut improves the signal to 
background ratio by approximately a factor of 1.5. 

Finally, for the HighTower-triggered data the requirement is made that at 
least one of the reconstructed decay photons alone satisfies this trigger. This 
requirement is made to guarantee that the trigger efficiency is the same in both 
real and simulated data, as was already mentioned in section [4.5.31 
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Figure 5.4: The energy asymmetry of the photon pairs reconstructed in p + p data. 



5.3 Invariant mass distribution 

After cuts, the pairs of BEMC points are turned into 4-vectors by assuming 
that the decay photons originate from the reconstructed main vertex. For each 
TT*^ candidate, the pseudorapidity ry, the azimuth (j), the transverse momentum 
Pt, and the invariant mass M^-y (Eg. 15. 1( ) are calculated. In Figure [531 we show 
the 77, (j), Pt, and Mj^ distributions of the 7r° candidates in the p + p dataset. 
For the d + Au data these distributions look similar as those shown for p + p. 

The r] distribution shows the decrease of the calorimeter acceptance at the 
edges, because there it is likely than one of the decay photons escapes the 
calorimeter. The asymmetry is due to the fact that the calorimeter half-barrel 
is positioned asymetrically with respect to the interaction point. The structure 
seen in the (j) distribution reflects the azimuthal dependence of the calorimeter 
acceptance caused by failing SMD modules. 

In Figure I5.5f c) is shown the pr distribution of the photon pairs separately 
for the MinBias and HighTower datasets. It is seen that the HighTower triggers 



56 



CHAPTER 5. NEUTRAL MESON RECONSTRUCTION 




significantly increase the rate of pion candidates at large px- The pT-integrated 
invariant mass distribution in Figure I5.5f d) clearly shows the ir^ and rj peaks 
superimposed on a broad background distribution. This background has a com- 
binatorial and a low mass component. In the next two sections we will discuss 
each background component in detail. 



5.4 Combinatorial background 

The combinatorial background in the invariant mass distribution originates from 
pairs of photon clusters that are not produced in a single tt^ decay. To describe 
the shape of the combinatorial background, we use the event mixing technique, 
where photon clusters from two different events are combined. To mix only sim- 
ilar event topologies, the data were subdivided into the mixing classes based on 
the vertex position, BEMC multiplicity, and trigger type (MinBias, HighTower- 



5.4. COMBINATORIAL BACKGROUND 



57 




MiiiBias 

HighTower-1 

HighTower-2 



-%0 -40 -20 20 40 60 ^0 5 10 15 20 25 30 35 

z^.,,,.;,,^ [cm] Number of BEMC points 

Figure 5.6: Distributions of the number of events, as a function of vertex z (left) and number 
of BEMC points (right), obtained from the p + p dataset. In both figures are shown the 
intervals used to define the event mixing classes. 



1, and HighTower-2). In Figure [5?6l we show the p + p vertex and multiphcity 
distributions, and the bins defining the mixing classes. 

Figure 15.71 shows an example of an invariant mass distribution in the 
4 < pt < 5 GeV/c bin, obtained from the HighTower-1 p + p data, together 
with the combinatorial background obtained from the event mixing. The mixed 
event background distribution is normalized to the same-event distribution in 
the invariant mass region 0.9 < M^^ < 1.2 GeV/c^. In the bottom panel of this 
figure the background subtracted distribution is shown. 

It can be seen that there is still some residual background in the interval 
0.2 < My^ < 0.4GeV/c^, which could be caused by the fact that the mixing 
procedure does not fully take into account the correlation structure of the event. 
For example, an important source of particle correlations is the jet structure, 
which is not present in the sample of mixed events. In order to preserve jet- 
induced correlations, the jet axes in both events are aligned before mixing, as 
described below. 

To determine the {rj, (f>) position of the most energetic jet in every event, the 
standard STAR jet finding algorithm [70] was used. The mixed pion candidates 
were constructed by taking two photons from different events, where one of the 
events was displaced in rj and (f) by Ar/ = 772 ~ ?7i and Acj) = 4>2 — 4>i, respectively. 
Here 771,2 and 0i,2 are the jet orientations in the two events. 

In Figure \bM we show a schematic view of two superimposed events, where 
the jet axes are aligned. In order to minimize acceptance distortions, the events 
were divided into mixing classes in the jet rj coordinate. By mixing only events in 
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Figure 5.7: The same-event invariant mass distribution and the combinatorial background 
obtained from the random event mixing (top) and the background subtracted distribution 
(bottom). The shaded area in the top plot indicates the region where the mixed event back- 
ground is normalized to the data. 
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Figure 5.8: A schematic view of two superimposed events, wiiere the jet axes are ahgned. 



the same class, the shift Arj was kept smaller than 0.1. Because the calorimeter 
has a cylindrical shape, the shift in (p does not induce any significant acceptance 
distortion. 

However, a side effect of this procedure is that correlations are induced if 
there is no real jet structure, because the jet finding algorithm will then simply 
pick the most energetic track in the event. To reduce possible bias introduced 
by such correlations, we assume that a jet structure is associated with large pr 
pions but not with low pT pious. The combinatorial background is then taken 
as a pT-dependent linear combination of the distributions obtained by random 
mixing and jet-aligned mixing, 

B{M,pt) = A{pt)Bj{M,pt) + (1 - A{pt))Br{M,pt). 

Here Bj{M,pt) and Bii{M,pT) are the background spectra from, respectively, 
the jet-aligned and random event mixing in a given px bin. The interpolation 
coefficient A{pt) is given by 

A{pT) = / f ^ + ^ fo//^ < 10 GeV/c 
^ ' \ 1 otherwise, ^ ' 

where the coefficients are a = 0.097 (GeV/c)~^ and h = —0.117. We assign a 
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systematic uncertainty of 10% to A, which propagates into a systematic uncer- 
tainty of 5% on the tt" and 1% on the rj yields. 

In Figure 15.91 we plot the same invariant mass spectrum as that shown in 
Figure 15. 7( with the background estimated by the combined random and jet- 
aUgned event mixing. The mixed event background is normalized to same-event 
distribution in the ranges 0.3 < M^^ < 0.4 and 0.8 < M^^ < 1.6GeV/c. By 
changing the subtracted background within the normalization uncertainty we 
obtained a systematic error on the tt^ and rj yields. This error was found to 
increase with px from 0.5 to 3% for the tt'^ and from 10 to 50% for the t] yield. 

In the bottom panel of Figure 15.91 the background subtracted spectrum is 
plotted, which still shows a residual background component at low invariant 
mass. The origin of this background is described in the next section. 

5.5 Low- mass background 

In Figure 15.11 we have shown a double peaked hit pattern, which will be recon- 
structed by the clustering algorithm as two separate adjacent clusters. How- 
ever, it is possible that random fluctuations will accidentally generate such a 
two peak structure, so that the clustering algorithm will incorrectly split the 
cluster. These random fluctuations enhance the yield of pairs with minimal 
angular separation and thus contribute to the lowest di-photon invariant mass 
region, as can be seen in Figure [5T9l However, at a given small opening angle the 
invariant mass increases with increasing energy of the photons, so that the low 
mass background spectrum will extend to larger values of Mjj with increasing 
Pt of the parent particle. 

The shape of the low mass background was obtained from a simulation as fol- 
lows. Single photons were generated with flat distributions in (f>, —0.2 < 77 < 1.2 
and < Pt < 25 GeV/c. These photons were tracked through a detailed 
description of the STAR geometry with the GEANT program [7T] . A detailed 
simulation of the electromagnetic shower development in the calorimeter was 
used to generate realistic signals in the towers and the SMD. The simulated 
signals were processed by the same reconstruction chain as the real data. 
Photons with more than one reconstructed cluster were observed, and the 
invariant mass and px of such cluster pairs were calculated. The invariant 
masses were histogrammed with each entry weighted by the pt spectrum of 
photons in the real data, corrected for the photon detection efficiency. 

In the top plot of Figure lS.lOl we show the low mass background distributions 
in three bins of the reconstructed pair px- It is seen that the distributions 
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Figure 5.9: The same-event invariant mass distribution and the combinatorial background 
obtained from the jet-aligned event mixing (top) and the background subtracted distribution 
(bottom). The shaded regions in the top plot indicate where the mixed event background is 
normalized to the data. 
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Figure 5.10: The low mass background distributions from erroneous splitting of single photons 
in three bins of the reconstructed pair (top). The distributions extend to larger invari- 
ant masses with increasing py. The low mass background component in the invariant mass 
distribution obtained from p + p data (bottom) . 
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indeed move to larger invariant masses with increasing px and extend far into 
the pion window at large pr- For this reason, it is not possible to estimate this 
background from a phenomenological fit to the data, so that we have to rely on 
the Monte Carlo simulation to subtract the low mass background. 

The second significant source of neutral clusters in the calorimeter are the 
neutral hadrons produced in the collision, mostly antineutrons. As a first 
attempt to account for the additional low mass background from these hadrons, 
simulations of antineutrons were performed in the same way as photons, and the 
reconstructed invariant mass distribution was added according to the realistic 
proportion n/7 = n/27r°. The ratio u/tt^ was taken to be equal to the average 
value of p/tt^ from the STAR measurement |28) in the px range covered by 
each of MinBias and HighTower datasets. In the bottom plot of Figure [5^. 101 we 
compare the simulated low mass background (histogram) to the data. 

In Figure 15.111 we show the invariant mass spectra and the low mass 
background component (top), together with the final background subtracted 
spectrum (bottom). 

5.6 Yield extraction 

The complete set of invariant mass spectra for all px bins, triggers, and datasets 
are shown in Figures [5 . 1 2H5 . 1 7l For display purposes, the spectra are normalized 
to the bin content in the tt'^ peak. The shaded areas in the figures indicate the 
TT^ and rj peak regions where the yields are calculated simply by adding up the 
bin contents. 

The left border of the 7r° peak region was taken to be a linear function of px, 
common for all datasets and triggers. It was adjusted in a way that most of the 
yield is captured, while the low mass background and its associated uncertainty 
is avoided as much as possible. The right border also linearly increases with 
px, in order to cover the asymmetric right tail of the peak. Similarly, the ry 
peak region is a pT-dependent window that captures most of the signal. For 
completeness, we give below the parametrization of the 7r° and ij windows: 

75 + l.lpx < M^o < 250 + 3.3pT MeV/c^, 
350 + 3.3pT < < 750 MeV/c^. 

The stability of the yields was determined by varying the vertex position cut, 
the energy asymmetry cut, and the yield integration window. From the observed 
variations, a point-to-point systematic error of 5% was assigned to the yields. 



64 



CHAPTER 5. NEUTRAL MESON RECONSTRUCTION 



g 



1200 
1000 
800 
600 
400 
200 

a 



- ' t 



— + 



S 1200 
o 

" 1000 
800 
600 
400 
200 

a 



p—p HighTowcr-1 
4 < 5 GoV/c 



t * 

+ + 



0.2 



0.4 



0.6 



0.8 



M^^ [GeV/c2 



Figure 5.11: The invariant mass distribution before (top) and after the low mass background 
subtraction (bottom). 




Figure 5.12; Invariant mass distributions in all pT bins, p + p MinBias data. The spectra are 
normalized to the bin content in the tt" peak. 



66 



CHAPTER 5. NEUTRAL MESON RECONSTRUCTION 




Figure 5.13; Invariant mass distributions in all pT bins, p + p HighTower-1 data. The spectra 
are normalized to the bin content in the peak. 
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Figure 5.14: Invariant mass distributions in all PT bins, p + p HighTower-2 data. The spectra 
are normalized to the bin content in the tt" peak. 




Figure 5.16; Invariant mass distributions in all px bins, d+An HigliTower-l data. The spectra 
are normalized to the bin content in the vr" peak. 
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Figure 5.17; Invariant mass distributions in all pT bins, d+ Au HighTower-2 data. The spectra 
are normalized to the bin content in the vr" peak. 



Chapter 6 

Invariant yield calculation 



The invariant yield of the neutral pions and r] mesons per one minimum bias 
collision, as a function of the transverse momentum pT, is given by 

^ d^N _ d^N _ d'^N 
dp^ Prdprdydcl) 2TTpTdpTdy' 

where in the last equation isotropic production in azimuth is assumed. Using 
the experimentally measured quantities, the invariant yield is calculated as 

^d^^^ gycncx AY 1 1 1_ 

dp3 27rpT iVtrigi^trig(l " £beam) ApT At/ £acc ^cpv ^j-y/^ ' 

where: 

• AY is the raw yield measured in the bin Apr At/; 

• -^trig is the number of triggers recorded; 

• Ktiig is the trigger prcscalc factor that is unity for the MinBias events and 
larger than unity for the HighTower data. The product A^trig-K^trig then 
gives the equivalent number of minimum bias events that produced the 
yield AY; 

• ^vertex is the vertex finding eSiciency in minimum bias events; 

• fi^beam is the beam background contamination in minimum bias events; 

• ApT is the width of the pt bin for which the yield is calculated; 
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• Ay is the rapidity range of the measurements, in this analysis Ay = 1; 

• £acc is the BEMC acceptance and efficiency correction factor; 

• £cpv is a correction for random vetoes; 

• F-y^/r is the branching ratio of the di-photon decay channel, equal to 0.988 
for ttO and 0.392 for 77 [55]. 

Each of these corrections is described in detail in one the following sections. 

6.1 Acceptance and efficiency correction 

To calculate the acceptance and efficiency correction factor eacc, a Monte Carlo 
simulation of the detector was used, where neutral pions and their decay photons 
were tracked through the STAR detector geometry using GEANT [7T]. The 
simulated signals were passed through the same analysis chain as the real data. 

The pions were generated in the pseudorapidity region --0.3<r7< +1.3, 
which is sufficiently large to account for edge effects caused by the calorimeter 
acceptance limits of < 77 < 1. The azimuth was generated flat in — tt <(/)<+ tt. 
The pt distribution was taken to be flat between zero and 25 GeV/c, which 
amply covers the measured pion px range of up to 17 GeV/c. The vertex 
distribution of the generated pions was taken to be Gaussian in z, with a spread 
of cr = 60 cm and centered at z = 0. 

The generated pions were allowed to decay into 7r° 77. The GEANT 
simulation accounts for all interaction of the decay photons with the detector, 
such as pair conversion into e+e^ and showering in the calorimeter or in the 
material in front. 

To reproduce a realistic energy resolution of the calorimeter, an additional 
smearing has to be applied to the energy deposit generated by GEANT in the 
towers. The effect of this can be seen in Figure 16.11 where the simulated 7r° 
invariant mass peak is shown in comparison to the p + p data with and without 
smearing. An additional spread of 5% was used to reproduce the p + p data and 
10% for the d + Au data. 

To reproduce the px spectrum of pions in the data, each Monte Carlo event 
was weighted by a pT-dependent function. Such weighting technique allows to 
sample the whole pt range with good statistical power, while, at the same time, 
the bin migration effect caused by the finite detector energy resolution is repro- 
duced. A next-to-leading order QCD calculation [72] provided the initial weight 
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Figure 6.1: The invariant mass distribution in the real p+p data (crosses) and in the simulation 
(histogram). The Monte Carlo produces a narrower tt" peak (left) than is observed in the data, 
so that an additional energy smearing was introduced to reproduce the calorimeter resolution 
(right). 



function, parametrized as described in Section 17.11 which was subsequently 
adjusted in an iterative procedure. 

As mentioned in Section 14. 3[ the time dependence of the calorimeter 
acceptance is stored in data tables, which are fed into the analysis. In or- 
der to reproduce this time dependence in the Monte Carlo, the simulated events 
were assigned time stamps that follow the timeline of the real data taking. In 
Figure 16.21 is shown, separately for MinBias and HighTower data, the accumu- 
lated real data statistics per day (histogram) , together with the time distribution 
of the simulated events (full circles). In this way, the geometrical calorimeter 
acceptance (fraction of good towers) was reproduced in the Monte Carlo with a 
precision of better than 0.5%. 

In the real data analysis, we use vertices reconstructed from the TPC tracks 
with a sub-millimiter resolution, as well as vertices derived from the BBC time 
of flight measurement with a precision of about 40 cm. To account for this poor 
resolution, a fraction of the simulated pions had their point of origin artificially 
smeared in the z direction. This fraction was taken to be 35% of the generated 
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Figure 6.2: Statistics accumulated per day (histogram) and simulated in the Monte Carlo (full 
circles), for p + p MinBias (left) and HighTower data (right). 



pions in case of the p+p MinBias analysis, and taken to be zero for all the other 
datasets since no BBC vertex was used in these sets (see Chapter 

In Figure [^751 we show the 77 and distributions of the reconstructed Monte 
Carlo pions in comparison to the p + p data. The agreement is satisfactory, 
indicating that the calorimeter acceptance is well reproduced in the simulation. 
In Figure 16.41 the reconstructed pT of simulated pions is compared to that of 
pion candidates from the p + p data. It is seen that the HighTower trigger 
threshold effects are reasonably well reproduced. 

In Figure IG.Sf a) the background subtracted invariant mass distribution is 
shown in the region 4 < < 5 GeV/c obtained from the p + p HighTower-l 
data, together with the corresponding distribution from the Monte Carlo. In 
order to compare the real and simulated invariant mass distributions for all 
bins in pT and for all datasets, we have estimated the position and width of 
the peaks by Gaussian fits in the peak region. In Figure IG.Sf b) are shown the 
peak positions obtained from the fit to the p + p data. It is seen that the 
peak position shifts towards higher masses with increasing pT- This shift is a 
manifestation of bin migration effects that originate from statistical fluctuations 
in the calorimeter response. Due to the steeply falling pT spectrum the energy 
resolution will cause a net migration towards larger pT- Since larger values of 
Pt imply larger values of M-y^, the migration effect will bias the invariant mass 
peak towards larger values. The good agreement between the data and Monte 
Carlo indicates that such resolution and migration effects are well reproduced. 

In Figure I6.5f c) is shown the comparison of the 7r° peak width in data and 
simulation. The peak width is well reproduced in simulation, which is not 
surprising since additional smearing was introduced to improve the comparison 
between data and Monte Carlo, see Figure [STTl 
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Figure 6.3: Distributions of rj (top) and <p (bottom) coordinates of the reconstructed Monte 
Carlo pions, compared to the pion candidates in the p + p HighTower-1 data. The structure 
seen in the </> distribution reflects the azimuthal dependence of the calorimeter acceptance, 
caused by failing SMD modules. This structure is well reproduced in the simulations. 
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[GeV/c] 

Figure 6.4: Distributions of the reconstructed of the Monte Carlo MinBias, HighTower-1, 
and HighTower-2 pions, compared to the p + p data. 



The acceptance and efficiency correction factor was calculated from the 
Monte Carlo simulation as the ratio of the raw yield of neutral pions recon- 
structed in a pt bin, to the number of simulated pions with the true pt in that 
bin. This was done separately for each trigger, using the same pion reconstruc- 
tion cuts as was done in the real data analysis. In particular, the reconstructed 
value of pseudorapidity was required to fall in the range < 77 < 1 in both the 
data and the Monte Carlo, while in the latter the generated value of 77 was also 
required to fall in this range. 

In Figure 16.61 arc shown the 7r° and ij correction factors for all datascts and 
triggers used in this analysis. 

The large difference between the MinBias and HighTower correction factors 
is caused by the SMD requirement in the HighTower data, while in the MinBias 
data we accept all reconstructed BEMC points. The absense of the SMD 
information also reduces the 7r° reconstruction efficiency at pt > 3 GeV/c, where 
the decay photons are separated by less than two towers. The t] reconstruction 
starts being affected at larger values of px- 
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Figure 6.5: Invariant mass spectrum reconstructed in the simulation, in comparison to the 
p + p HighTower-1 data in 4 < pT < 5 GeV/c bin (top). Peak position (middle) and width 
(bottom) in the real data and MC simulation. 
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Figure 6.6: Acceptance and efficiency factor Eacc calculated from the Monte Carlo simulation 
for the p+p (left-hand plots) and d+ Au datasets (right-hand plots). The tt" and r} efficiencies 
are shown separately in the top and bottom plots, respectively. 



0.14 
0.12 
0.1 
0.08 
0.06 
0.04 
0.02 







p—p HighTower- 1 vr" 
Relaxed SMD cut 



10 



12 



14 16 
[GcV/c] 



Figure 6.7: Acceptance and efficiency correction for the p+p HighTower-1 data, with standard 
set of cuts and with SMD quality cut removed. 
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We have checked the effect of the SMD quality requirement (at least two 
adjacent strips in a cluster) on the correction factor for HighTower triggered 
data. In Figure 16.71 is shown the correction factor calculated for the p + p 
HighTower-1 dataset with (full squares) and without (crosses) the SMD quality 
requirement. It is seen that this requirement reduces the number of accepted tt'^ 
candidates by about 45%. This explains the difference between the HighTower-1 
and IIighTower-2 (no SMD quality cut) correction factors at large px seen in 
Figure 16.61 

To verify a possible dependence of the acceptance correction on the track 
multiplicity and thus on the centrality, we have analyzed a sample of generated 
neutral pions embedded in real d + Au data. These embedded data are cen- 
trally produced by the STAR offline group and are used by several analyses in 
STAR [73]. No significant centrality dependence was found, so that same cor- 
rection factors were applied to the various centrality classes in the d + Au data. 

6.2 Corrections for random vetoes 

This analysis uses the TPC as a veto detector to reject charged particles, which 
introduces false rejection of photon clusters if an unrelated charged particle 
happens to hit the calorimeter nearby the cluster. In Figure 16.81 we plot the 
distribution of distances between the BEMC point and the closest charged track 
in the event. In this plot one easily distinguishes the peak of real charged 
particles at small distances, superimposed on a random component, which shows 
up as a shoulder at larger distances. Assuming that the charged tracks are 
uniformly distributed in rj and (j) around the BEMC point, it follows that the 
radial distribution is given by 

fiD)^De-^^, (6.3) 

where the parameter t has the meaning of the local track density in the region 
where the photon probes it. This parameter is obtained from a simultaneous fit 
to the data in all bins of the event multiplicity M, assuming its linear dependence 
on the multiplicity t ^ a + bAI. The parametrization (|6.3p well describes the 
random component, as shown by the full curve in Figure 16.81 The relative 
amount of random coincidences is then obtained by integrating the fitted curve 
up to the distance cut and weighting with the multiplicity distribution observed 
in each px bin. Separate sets of correction factors were calculated for the various 
triggers in the p + p and d + Au data. The results are shown in Figure 16.91 as 
a function of px. We have applied a correction factor of £cpv — 0.94 ± 0.02 to 
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D =\lArf I A<f)'^ 

Figure 6.8: Distribution of the distance between BEMC points and the closest track, obtained 
from p + p HighTower-1 data in the bin 4 < pj- < 5 GeV/c. The curve shows a fit to Eq. (16.31 1 
and the vertical line indicates the CPV cut. 



the p+p datasets and of Scpv = 0.89 ± 0.02 to the d + Au datasets. The errors 
assigned to these corrections contribute to a pt independent systematic error 
on the corrected tt" and 77 yields. 



6.3 High Tower trigger scale factors 

We have shown in Figure IS.Sf c) the pr distribution of tt'^ candidates for the 
p+p MinBias, HighTower-1, and HighTower-2 data. To match the HighTower 
spectra to those of the MinBias, a pT-independent scale factor (i^trig) was ap- 
plied. These scale factors were estimated as the ratio of observed MinBias to 
HighTower event rates 
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Figure 6.9: Charged particle veto correction in p + p (left) and d + Au (right) data. The 
horizontal line indicates the correction factor £cpv applied to the data, while the shaded band 
corresponds to the systematic uncertainty assigned to the correction factors. 



Here A^mb and A^ht are the numbers of MinBias and HighTower triggers that 
pass the event selection cuts described in Chapter S) The factors 5mb and 5ht 
are the hardware prescale factors adjusted on a run-by-run basis to accomodate 
the DAQ bandwidth. In Eq. ()6.4p . the sums are taken over aU runs in which 
both the MinBias and HighTower triggers were active. 

To check the resuhs, the scale factors were also estimated using another 
method. Here the HighTower software filter (see section 14.5.31) was applied to 
the minimum bias data. The scale factors were then obtained as the ratio of the 
total number of MinBias events to the number of those that passed the filter. 
To obtain a more precise HighTower-1 / HighTower-2 relative scale factor, the 
software filter was applied to the HighTower-1 dataset. 

The results from the two methods agree within 3% for HighTower-1 data and 
within 5% for HighTower-2 data. This is taken as the systematic uncertainties 
on the trigger scale factors. 



6.4 Vertex finding efficiency 

In the tt'^ reconstruction it is assumed that the decay photons originate from the 
vertex. It is therefore required that each event entering the analysis has a recon- 
structed vertex. In the p+p dataset this requirement is always fulfilled, because 
we use the BBC timing information in case the TPC vertex reconstruction fails 
(this happens in about 35% of the minimum bias events). 

In the d + Au HighTower data, the charged track multiplicities are large 
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enough to always have a reconstructed TPC vertex. However, a TPC vertex 
is missing in about 7% of the minimum bias events and cannot be recovered 
from BBC information because the BBC is not included in the d + Au mini- 
mum bias trigger. Minimum bias events without vertex have low charged track 
multiplicity, and the contribution from these very soft events to the tt" yield 
above 1 GeV is assumed to be negligible [71]. Therefore, the correction for ver- 
tex inefficiency is applied as a constant normalization factor to the yield and its 
uncertainty contributes to the total normalization uncertainty of the measured 
cross sections. 

The vertex efficiencies were determined to be 0.93 ± 0.01 from a full simula- 
tion of d + Au minimum bias events as described in |74| . However, this efficiency 
depends on the centrality, and we assume that central events are 100% efficient. 
Scaling the above efficiency by the ratio of peripheral to total number of d + Au 
events, we obtain an efficiency correction factor of 0.88 ± 0.02 for the sample of 
peripheral events. 

Note, that the difference between vertex finding efficiencies in MinBias and 
HighTower data is effectively absorbed in the scale factor i^trig defined in the 
previous section. The vertex finding efficiency correction is, therefore, applied 
to the minimum bias data, as well as to the scaled HighTower trigger data. 



6.5 Residual beam background contamination 

The beam background contamination in the d + Au minimum bias trigger has 
been estimated from an analysis of the RHIC empty bunches to be 5 ± 1 % [7S] . 
In our analysis, the beam background in d + Au events is rejected when the 
energy deposit in the calorimeter is much larger than the total energy of all 
charged tracks reconstructed in the TPC, see Section 14.51 To estimate the 
residual beam background in our data, we have analysed a sample of 3 x 10^ 
minimum bias triggers from unpaired RHIC bunches. These events were passed 
through the same analysis cuts and reconstruction procedure as the real data. 
We observed that about 10% of the fake triggers passed all cuts, and that 
none of these contained a reconstructed tt'^. The residual beam background 
contamination is thus estimated to be 0.1 x 5 = 0.5%, which is considered to be 
negligible. 

In the p + p data the beam background contamination to the minimum 
bias trigger rate is also estimated to be negligible due to the BBC coincidence 
requirement in the trigger and the cut on the BBC vertex position. 
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6.6 Bin centering scale factors 

To assign a value of pT to the yield measured in a pT bin, the procedure from |76| 
was applied. Here the measured yield, initially plotted at the bin centers, is 
approximated by a power law function of the form 

fiPr) = 77- i—rw- (6.5) 

{1+Pt/po) 

To each bin a momentum p^ was assigned as calculated from the equation 

/(Pt) = / /(^)da;- (6-6) 

^PT J Apt 

The function (j6.5|) is then re- fitted taking p^ as the abscissa. This procedure was 
re-iterated until the values of p^ were stable (typically after three iterations). 
Final fitted curves are shown in Figures 16.111 and 16.121 

For convenience of comparing results from the various datasets, the yields 
were scaled to the px bin centers by the ratio 

K = f6 7) 

where px is the center of the bin. The statistical and systematic errors were 
also scaled by the same factor. 



6.7 Jacob ian correction 

All calculations in this analysis were performed in the defined pesudorapidity 
region < 77 < 1 that corresponds to the rapidity region < y < yo, where 
the rapidity limit yo is well approximated by pseudorapidity for a particle with 
momentum much larger than its mass. 

The correction was applied to account for the rapidity limit yo being not 
equal to pseudorapidity 77 = 1, as shown in Figure 16.101 This correction is 
smaller than 10% for the 77 data points at pT < 3 GeV/c, and is negligible for 
the other data points. 



6.8 Fully corrected yields 

The fully corrected tt*^ invariant yields per minimum bias event in p + p and 
d + Au collisions were calculated from Eq. (|6.2p and are shown in the top plots 
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0.95 




[GeV/c] 

Figure 6.10: Jacobian correction that accounts for the rapidity hmit yo being not equal to 
pseudorapidity rj = I. 



of Figure 16.111 and Figure 16.121 The curves in these figures represent a fit of 
Eq. (|6.5p to the data. In the bottom plots are shown the ratios between the 
data and the fit. From these plots it is seen that the agreement between the 
datasets taken with the three triggers is satisfactory. 

For the calculation of the final cross section results and cross section ratios, 
the data from three triggers were merged together and only one data point was 
chosen in each overlapping pT bin. The HighTower-1 points were preferred over 
MinBias, and HighTower-2 over HighTower-1, because at high pT data samples 
are highly correlated, while HighTower datasets typically have smaller statistical 
error. 

The systematic uncertainty due to the calorimeter calibration was estimated 
from 

df 



SfiPT) 



dpq 



Spi 



where Spx was taken to be 5% in the d + Au and p + p data (see Section 
and where the derivative was calculated from the fitted function, Eq. ()6.5|) . This 
PT-dependent systematic uncertainty is, on average, 38% in the p + p data and 
44% in the d + Au data. 
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Figure 6.11: Invariant yield of tt" per minimum bias event in p + p and d-l- Au collisions (top). 
Curves are the power law fits given in the text. Invariant yield divided by the fit to the p -f p 
(middle) and d + Au (bottom) data. The errors shown are statistical only. 
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Figure 6.12: Invariant yield of ij meson per minimum bias event in p + p and d + An collisions 
(top). Invariant yield divided by the fit to the p + p (middle) and (i + Au (bottom) data. The 
errors shown are statistical only. 
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All systematic error contributions mentioned in this and in the previous 
sections are summarized in Table [01 classified into the following categories: 

• A point-by-point systematic uncertainty; 

• B point-by-point p^-correlated systematic uncertainty, but uncorrelated 
between datasets; 

• C point-by-point pr-correlated systematic uncertainty, also correlated 
between datasets; 

• N normalization uncertainty, uncorrelated between datasets. 



Table 6.1: Systematic error contributions. The classifications A, B, C and N 
are defined in the text. The error contributions to the cross section, the rj/ir'^ 
ratio, Rcp and RdA are indicated in the respective columns. The last column 
refers to the section where each source of systematic error is described. 

Source Type E d^cr/dp^ ri/n^ Rcp RdA Section 



Combinatorial background 


A 


+ 


+ + 


h El 


Mixed-event background 


C 


-t- 




[531 


Random vetoes 


N 




+ 


h [O 


HighTower scale factors 


B 


+ 


+ 


h E3] 


Analysis cuts 


A 


+ 


+ + 


h ESI 


Energy scale 


B 


+ 




h 


Vertex finding efficiency 


N 




+ 


h m 


Min. bias cross section 


N 


+ 




h o 


Glauber model A'coU 


N 




+ 
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Chapter 7 

Results and discussion 



7.1 Cross section 



The invariant differential cross section for 7r° and ij production in inelastic p + p 
interactions is given by 

dp3 dp3 -''NSD2^p^^p^^y- l^-iJ 

It has been shown that the singly diffractive (SD) contribution to the inelastic 
cross section is negligible at > 1 GeV/c [77| , so that we can assume that the 
differential inelastic cross section is equal to the differential NSD cross section 
in our pT range. The total NSD cross section in p + p collisions was taken to 
be cr^SD = ^'-'•0 ^ "^-^ ^ described in Section HTT] The total hadronic cross 
section in d + Au collisions was taken to be cr^+^" = 2.21 ± 0.09 b [H]. 

In tables 17.11 and 17.21 we list the cross sections calculated from Eq. (|7.ip 
for the p + p and d + Au datasets. In the third column of these tables are 
given the statistical errors, while in the remaining columns the quadratic sum 
of the systematic errors are given separately for each group defined in Table EUl 
In addition to these pr-dependent systematic errors, the quadratic sum of the 
normalization uncertainties is found to be 12.2% for the p + p and 5.6% for the 
d + Au data. 
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Table 7.1; Invariant cross section of tt^ production measured inp + p collisions. 
Systematic errors classification given in Section [6.81 Normalization uncertainty 
of 12.2% is not included. 

Pt E d?a/Ap^ Statistical Systematic errors 



[GeV/c] 


[mb GeV"^ c^] 


error 




A 




B 




C 




1.25 


2.646x10" 


-1 


5.973x10" 


-3 


1.336x10" 


-2 


9.835x10" 


-2 


1.323x10" 


-2 


1.75 


5.095x10" 


-2 


1.217x10" 


-3 


2.579x10" 


-3 


1.898x10" 


-2 


2.548x10" 


-3 


2.25 


1.314x10" 


-2 


4.349x10" 


-4 


6.669x10" 


-4 


4.902x10" 


-3 


6.569x10" 


-4 


2.75 


4.154x10" 


-3 


1.948x10" 


-4 


2.115x10" 


-4 


1.553x10" 


-3 


2.077x10" 


-4 


3.50 


5.724x10" 


-4 


6.062x10" 


-5 


2.928x10" 


-5 


2.146x10" 


-4 


2.862x10" 


-5 


4.50 


6.076x10" 


-5 


2.717x10" 


-6 


3.131x10" 


-6 


2.293x10" 


-5 


3.038x10" 


-6 


5.50 


1.223x10" 


-5 


4.403x10" 


-7 


6.353x10" 


-7 


4.632x10" 


-6 


6.113x10" 


-7 


6.50 


3.246x10" 


-6 


1.248x10" 


-7 


1.702x10" 


-7 


1.234x10" 


-6 


1.623x10" 


-7 


7.50 


9.592x10" 


-7 


2.638x10" 


-8 


5.081x10" 


-8 


3.681x10" 


-7 


4.796x10" 


-8 


9.00 


2.362x10" 


-7 


6.919x10" 


-9 


1.272x10" 


-8 


9.114x10" 


-8 


1.181x10" 


-8 


11.00 


5.797x10" 


-8 


3.029x10" 


-9 


3.198x10" 


-9 


2.253x10" 


-8 


2.898x10" 


-9 


13.00 


1.632x10" 


-8 


1.564x10" 


-9 


9.250x10" 


-10 


6.389x10" 


-9 


8.162x10" 


-10 


15.50 


4.357x10" 


-9 


5.918x10" 


-10 


2.559x10" 


-10 


1.720x10" 


-9 


2.178x10" 


-10 



Table 7.2: Invariant cross section of tt*' production measured in rf+Au collisions. 
Systematic errors classification given in Section [6.81 Normalization uncertainty 
of 5.6% is not included. 



Pt E d?a/Ap^ Statistical Systematic errors 
[GeV/c] [mb GeV~^ e^] error A B C 



1.25 


8.487x10^ 




3.716x10° 




4.286x10° 




3.133x10^ 




4.244x10° 




1.75 


2.052x10^ 




9.744x10" 


-1 


1.039x10° 




7.648x10° 




1.026x10° 




2.25 


6.116x10° 




3.602x10" 


-1 


3.104x10" 


-1 


2.302x10° 




3.058x10" 


1 


2.75 


1.643x10° 




1.556x10" 


-1 


8.365x10" 


■2 


6.245x10" 


1 


8.215x10" 


2 


3.50 


2.709x10" 


1 


4.373x10" 


■2 


1.386x10" 


■2 


1.044x10" 


1 


1.354x10" 


2 


4.50 


3.000x10" 


■2 


3.898x10" 


3 


1.546x10" 


3 


1.182x10" 


2 


1.500x10" 


3 


5.50 


6.924x10" 


3 


5.282x10" 


-4 


3.598x10" 


■4 


2.778x10" 


3 


3.462x10" 


■4 


6.50 


1.573x10" 


3 


1.002x10" 


-4 


8.249x10" 


-5 


6.457x10" 


■4 


7.864x10" 


■5 


7.50 


4.717x10" 


■4 


2.663x10" 


-5 


2.499x10" 


■5 


1.971x10" 


4 


2.359x10" 


■5 


9.00 


1.014x10" 


■4 


5.453x10" 


-6 


5.462x10" 


■6 


4.347x10" 


5 


5.071x10" 


6 


11.00 


2.439x10" 


■5 


2.172x10" 


-6 


1.346x10" 


■6 


1.081x10" 


5 


1.220x10" 


6 


13.00 


7.896x10" 


■6 


1.118x10" 


-6 


4.475x10" 


■7 


3.614x10" 


■6 


3.948x10" 


■7 


15.50 


1.585x10" 


■6 


4.056x10" 


-7 


9.311x10" 


■8 


7.542x10" 


7 


7.926x10" 


■8 
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To parametrize the px dependence, the measured 7r° cross sections were 
fitted to the power law function from Eq. (|6.5p . resuhing in the following 
parameters: 

Dataset A po n 

[mb GeV-^ c^] [GeV/c] 

p + p 7.53x102 0.95 9.31 
d + AvL 4.10 xlO^ 1.60 10.43 

The measured cross sections for tt'^ production in p+p collisions are shown in 
Figure Fm compared to the NLO pQCD calculation from Rcf. TS]. Input to this 
calculation are the CTEQ6M parton densities |3lJ and the KKP fragmentation 
functions [33 . The factorization scale fi was set equal to px and was varied 
by a factor of two to estimate the scale uncertainty, as indicated by the shaded 
band in the bottom plot of Figure [771] that shows the ratio of the measured cross 
sections to the QCD prediction. 

The errors shown in the plot are the statistical and point-to-point systematic 
uncertainties added in quadrature, excluding the uncertainty due to the energy 
calibration of the calorimeter. This additional uncertainty is shown by the outer 
lines around the data points on the lower plot. The normalization uncertainty 
is indicated by shaded band around unity on the right hand side of the plot. 

The tt" cross section measured in d + Au collisions is shown in Figure 17.21 
and compared to the NLO pQCD calculations of [72 . Here were used the KKP 
fragmentation functions, the CTEQ6M parton distributions for deuterium and 
the nuclear parton distributions for Au [79l [80l [81] . The errors shown in the 
plot are defined in the same way as in Figure [7?T] for p + p. 

It is seen that the measured tt'^ cross section in both the p + p and d + Au 
collisions is well described by the pQCD calculations in the full px range. The 
possible excess relative to the theory seen in the d + Au data at low px may 
be an indication of the Cronin effect, which was not included in the pQCD 
calculations. 

In Figure 17.31 we compare the tt" measurement in p + p and d + Au with 
previous measurements of charged pious by STAR [5S] . For ease of comparison, 
the TT* data points are divided by the 7r° pQCD curves. Note, that the normal- 
ization uncertainty shown by the grey bands in the figure are largely correlated 
between the and the tt^ data points. It is seen that the neutral and charged 
pion spectra agree with each other very well in both p + p and d + Au datasets. 

In Figure[73|we compare the present tt^ measurements with the neutral pion 
results from PHENIX [751 [20]. This comparison indicates a good agreement, 
within errors, between the results of the two experiments. 
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Figure 7.1: Cross section of the neutral pion production in p + p collisions (top), divided 
by the pQCD calculation (bottom). The errors are statistical and point-to-point systematic, 
excluding the energy calibration uncertainty shown as the outer lines (bottom). Normalization 
uncertainty is indicated by a shaded band around unity (bottom). 
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Figure 7.2: Cross section of the neutral pion production in d + Au collisions (top), divided 
by the pQCD calculation (bottom). The errors are statistical and point-to-point systematic, 
excluding the energy calibration uncertainty shown as the outer lines (bottom) . Normalization 
uncertainty is indicated by a shaded band around unity (bottom). 
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Figure 7.3: Cross section of the neutral pion production in p + p (top) and d + Au (bottom) 
collisions divided by the pQCD calculation, compared to the STAR tt^ 1281 . The errors 
are statistical and point-to-point systematic, excluding the calorimeter energy calibration 
uncertainty shown as the outer lines. Common normalization uncertainty is indicated by a 
shaded band around unity. 
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Figure 7.4: Cross section of the neutral pion production in p + p (top) and d + Au (bottom) 
collisions divided by the pQCD calculation, compared to the PHENIX tt" measurements in 
p + p 1781 and d + Au :20 . The errors are statistical and point-to-point systematic, excluding 
the energy calibration uncertainty shown as the outer lines. Normalization uncertainty is 
indicated by a shaded band around unity. 
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Table 7.3: rj/n'^ ratio measured inp+p collisions. Systematic errors classification 
given in Section 



Pt 




Statistical 


Systematic error 


[GeV/c] 




error 


A 


C 


1.25 


0.308 


0.069 


0.045 


0.012 


1.75 


0.380 


0.075 


0.059 


0.015 


2.25 


0.503 


0.097 


0.084 


0.020 


2.75 


0.292 


0.107 


0.052 


0.012 


3.50 


0.545 


0.166 


0.107 


0.022 


4.50 


0.599 


0.054 


0.131 


0.024 


5.50 


0.477 


0.045 


0.116 


0.019 


6.50 


0.593 


0.064 


0.158 


0.024 


7.50 


0.378 


0.038 


0.110 


0.015 


9.00 


0.381 


0.044 


0.125 


0.015 


11.00 


0.263 


0.076 


0.099 


0.011 


13.00 


0.343 


0.147 


0.146 


0.014 


15.50 


0.285 


0.188 


0.139 


0.011 



7.2 Eta to pion ratio 



The ry measurement is presented here as the ratio of i] to tt*' invariant yields, 
which allows many systematic uncertainties to cancel, see Table 16.11 in Sec- 
tion [6i8l The rj/n'^ ratios measured in p + p and d + Au collisions are listed in 
Tables 17.31 and 17.41 and are shown in Figure 17.51 The error definitions in the 
tables and in the plot are the same as described above for the differential cross 
sections. The present measurement agrees very well with previous PHENIX 
results [H2], as shown by the open squares in the plot. The full curves in 
Figure 17.51 show the asymptotic ratio R^/„o ~ 0.5 consistent with the world 
rj/ir^ measurements. The constant fit to our data at px > 4GeV/c gives 
= 0.42 ± 0.05 and R'^+^^' ^ 0.37 ± 0.08. The dotted curves in Figure O 
show the prediction based on empirical m^-scaling observation )83j that the 
hadron production cross sections have the same shape as a function of the 
transverse mass of the produced particle uit = -^/m^ + Pt- It is seen that the 
data are consistent with such scaling behaviour. 
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Figure 7.5: 77/77" ratio measured in p + p (top) and d + Au (bottom) collisions, compared 
the PHENIX measurements | 82| . Errors are statistical and systematic combined. 
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Table 7.4: rj/n^ ratio measured in d + Au collisions. Systematic errors classification 
given in Section 



Pt 




Statistical 


Systematic 


error 


[GeV/c] 




error 


A 




C 


1.25 


0.540 


0.344 


0.078 


0, 


,022 


1.75 


1.054 


0.310 


0.164 


0, 


,042 


2.25 


0.423 


0.337 


0.071 


0, 


,017 


2.75 


1.130 


0.459 


0.202 


0, 


,045 


3.50 


0.406 


0.601 


0.079 


0, 


,016 


4.50 


0.519 


0.355 


0.114 





,021 


5.50 


0.641 


0.189 


0.156 


0, 


,026 


6.50 


0.350 


0.175 


0.093 


0, 


,014 


7.50 


0.393 


0.134 


0.114 


0, 


,016 


9.00 


0.323 


0.109 


0.106 


0, 


,013 


11.00 


0.212 


0.168 


0.080 


0, 


,008 


13.00 


0.442 


0.254 


0.188 


0, 


,018 



7.3 Nuclear modification factor 

We calculate the RdA ratio, defined by Eq. (I1.9|l and (I1.10|) . as 

'■^ (A^coii) d^aP+P/dpTdy' ^'■^> 

where the nucleon-nucleon inelastic cross section (J^^ is taken to be 42 mb 
and (A^coii) = 7.5 ± 0.4 is calculated from the Glauber model, as described in 
Section IMl 

The nuclear modification factors for 7r° and rj are listed in Tables [731 and FTBl 
and shown in Figure [721 Again, the definition of the errors is as given for the 
differential cross sections in Section [^?51 Also shown in Figure [7TB1 are the results 
of RdA for charged pions measured by STAR \28^. A good agreement between 
STAR neutral and charged pions is observed. 

In Figure 17.71 we compare the RdA ratio for tt'^ (top panel) and rj (bottom 
panel) to the corresponding PHENIX measurements [HI] . Our data agree 
reasonably well with PHENIX, except at pt < 8 GeV/c, where the present 
results seem to be systematically higher by about 30%. 

The Rcp ratio for 7r° is listed in Table [777[ and shown in Figure 17751 compared 
to the STAR charged pions [28^. It is seen that the agreement between the 
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Table 7.5: Nuclear modification factor RdA for 7r°. Systematic errors classification 
given in Section [6751 Normalization uncertainty of 14.5% is not included. 



Pt 


RdA 


Statistical 


Systematic 


error 


[GeV/c] 




error 


A 




B 


1.25 


0.817 


0.040 


0.058 


0, 


,462 


1.75 


1.025 


0.055 


0.073 


0, 


.580 


2.25 


1.185 


0.080 


0.085 


0, 


.670 


2.75 


1.007 


0.106 


0.072 


0, 


.570 


3.50 


1.205 


0.233 


0.087 


0, 


.681 


4.50 


1.257 


0.173 


0.091 


0, 


.713 


5.50 


1.442 


0.122 


0.105 


0, 


.818 


6.50 


1.512 


0.104 


0.111 


0, 


.862 


7.50 


1.252 


0.079 


0.093 


0, 


,714 


9.00 


1.093 


0.067 


0.082 


0, 


.623 


11.00 


1.071 


0.111 


0.083 


0, 


,611 


13.00 


1.232 


0.211 


0.098 


0, 


,702 


15.50 


0.926 


0.268 


0.076 


0, 


,528 



Table 7.6: Nuclear modification factor R^a for V- Systematic errors classification 
given in Section [6.81 Normalization uncertainty of 14.5% is not included. 
Pt RdA Statistical Systematic errors 



[GeV/c] 




error 


A 


B 


2.75 


3.896 


2.088 


0.840 


2.204 


3.50 


0.898 


1.345 


0.218 


0.508 


4.50 


1.088 


0.736 


0.303 


0.617 


5.50 


1.939 


0.578 


0.612 


1.100 


6.50 


1.053 


0.525 


0.371 


0.600 


7.50 


1.302 


0.456 


0.507 


0.742 


9.00 


0.927 


0.325 


0.413 


0.529 
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• tt'' (this analysis) 

* STAR (7r+-7r-)/2, < 0.5 




[GcV/c| 

Figure 7.6: RdA ratio for vr'' (top) and rj meson (bottom), compared to the STAR 1281 . 
Errors are statistical and point-to-point systematic, excluding the calorimeter energy calibra- 
tion uncertainty shown as the outer lines. Common normalization uncertainty is indicated by 
a shaded band around unity. 
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Figure 7.7: RdA ratio for vr" (top) and t] meson (bottom), compared to the PHENIX mea- 
surements 82 84 . Errors are statistical and point-to-point systematic, excluding the energy 
calibration uncertainty shown as the outer lines. Normalization uncertainty is indicated by a 
shaded band around unity. 
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Table 7.7: Nuclear modification factor Rcp for tt". Systematic errors classification 
given in Section [^751 Normalization uncertainty of 11.1% is not included. 
Pt Rcp Statistical Systematic errors 



GcV/c] 






error 




A 




B 


1.25 


1 


032 


0.094 





073 





000 


1.75 


1 


177 


0.111 





084 





000 


2.25 


1 


265 


0.164 





090 





000 


2.75 


1 


059 


0.238 





076 





000 


3.50 


1 


211 


0.610 
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000 


4.50 


1 


428 


0.236 
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Figure 7.8: Rcp ratio measured in d + Au collisions, compared to STAR charged pions |28) . 
Errors are statistical and point-to-point systematic. Common normalization uncertainty is 
indicated by a shaded band around unity. 



neutral and charged pion measurements in STAR is very good. The ratio stays 
constant at a value consistent with unity beyond pp — S GeV/c. The indication 
of a decrease from the charged pion data is not supported by this measurement. 



7.4. CONCLUSIONS AND OUTLOOK 
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7.4 Conclusions and outlook 

There is a good agreement between the tt*^ cross sections in p + p and d + Au 
cohisions and nuclear modification factors measured in the present analysis, 
and the charged pions previousely measured in STAR at pT < 10 GeV/c. This 
demonstrates a consistency between the charged and neutral pion results, in 
spite of very different analysis methods and detectors (BEMC versus TPC) 
used for the measurements. This analysis extends the pT range of identified 
hadron measurements in STAR up to 17 GeV/c. There is also a good agreement 
with the corresponding tt" cross sections measured by the PHENIX experiment 
and with those calculated in NLO pQCD. 

From the measurement of the nuclear modification factor RdA, no suppression 
of the tt" production is seen in the d + Au collisions compared to p + p colli- 
sions. This is in line with the observation made elsewhere that the large 
suppression seen in the central Au + Au collisions is due to the final state effects. 

This analysis presents the first rj meson measurement in STAR. The cross 
section, presented as an 77/77° ratio, is in agreement with the PHENIX measure- 
ment and with the mT-scaling assumption. 

There are several important and unique features in the present analysis. 
First, the technique of estimating the low invariant mass background using the 
single photon simulation allows to remove the ir'^ contamination at high py, 
where this type of background is indistinguishable from the signal. A possi- 
ble further improvement would require a better handle on the SMD response 
simulation. Second, the jet-aligned event mixing method very well reproduces 
the combinatorial background in the tt^ and 77 peak region, so that no residual 
background subtraction is necessary. This eliminates the systematic uncertainty 
usually related to a residual background parametrization. 

It is seen that the experimental uncertainties can be significantly reduced 
by improving the calorimeter energy calibration. Better measurement in the 
low pt region would also require improvements in the simulation of the SMD 
response. Furthermore, calorimeter-based measurements using data taken since 
2006 will benefit from the full \ri\ < 1 BEMC acceptance coverage. 

In summary, these p + p and d + Au results provide a baseline measurement 
for the future Au + Au measurements. These measurements are interesting to 
shed light on quark number scaling in particle production at intermediate px 
and to study the origin of suppression phenomena at large pt- 
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Appendix A 

BEMC electronics 
operation 

The tower phototubes are powered by Cockroft-Walton (CW) bases that are 
able to keep the high voltage up to a high precision. The bases are programmed 
through the serial line from a dedicated computer in the Control Room. The 
analog signals fr'om the phototubes are routed to the tower digitizer crates 
mounted on the outer side of the magnet. 

The tower digitizer crate contains five boards that take 32 analog PMT 
inputs each and digitize it to 12 bit on each RHIC bunch crossing, storing in 
the digital pipeline until a level-0 trigger arrives. The crate controller board 
then sends the data packets to the Tower Data Collector on the platform that 
feeds it to the DAQ. The crate controller is also responsible for the slow control 
communication . 

The STAR level-0 trigger uses the BEMC data in the form of trigger prim- 
itives calculated by the tower digitizer boards, instead of the full tower data. 
Two trigger primitives are calculated for each tower patch of 0.2 x 0.2 in rj x ip 
(4x4 towers) using pedestal subtracted tower ADC: 

• High Tower - single largest tower signal in a patch 

• Patch Sum - sum of all 16 towers signal in a patch 

In the process of calculating those primitives, the on-board FPGA algorithm 
performs the following operations, as illustrated on Figure lA. II 
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APPENDIX A. BEMC ELECTRONICS OPERATION 



1. Drop the last 2 bits of the tower ADC, which becomes a 10 bit signal. 

2. Subtract the stored pedestal PED from the ADC, mask the channel out 
if necessary. The pedestals are calculated in a special way as described 
below. 

3. For the High Tower: convert 10 bits into 6, using one of four methods 
(HT6 selector), then select the largest value of all 16 towers as output. 

4. For the Patch Sum: drop the last 2 bits to make it an 8 bit value, then sum 
those from all 16 channels into a 12 bit value and transform it into the 6 
bit output value. The transformation function has a special shape which 
is described in details later. Internally, it uses a lookup table (LUT); the 
6 bit number stored in the LUT is the output. The PatchSum trigger 
sensitivity is, therefore, fixed: 

ADCP^*^'" = 16 X ADCtower. 



Finally, two 6 bit numbers are sent to the trigger Data Storage and Manipulation 
(DSM) boards upon recceiving a trigger signal. 

The following HT6 methods are available to select 6 of 10 bits for the High 
Tower output with varying degree of sensitivity, almost equivalent to selecting 
a constant attenuation factor: 

(0) Select 6 lowest bits, combine five highest bits by logical "and" into the 
highest bit of the result. This is the most sensitive trigger setting, one 
HighTower trigger ADC count is equal to 4 raw tower ADC counts. 

(1) Select 6 lower bits starting from 1, combine four highest bits by logical 
"and" into the highest bit of the result: ADC^g^"™' = 8 x ADCtower 

(2) Select 6 lower bits starting from 2, combine three highest bits by logical 
"and" into the highest bit of the result: ADCf/jlg^""^^' = 16 x ADCtower 

(3) Select 6 lower bits starting from 3, combine two highest bits by logical 
"and" into the highest bit of the result. This is the least sensitive trigger 
setting, one HighTower trigger ADC count is equal to 32 raw tower ADC 
counts. 

The tower pedestals and masks, the HT6 selectors, and the LUT arrays are 
prepared and uploaded into the on-board registers via the slow control program. 
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Figure A.l: The digital processing in the tower digitizer boards. 
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APPENDIX A. BEMC ELECTRONICS OPERATION 



The tower pedestals arc being specially prepared in a way that puts the 6 bit 
High Tower and Patch Sum pedestals at 1 (not zero) to be observed during the 
run. For each tower, the calculation starts from the exact value of the pedestal, 
which is measured by issuing the software triggers to FEE via slow control in a 
periods between data taking when there is no beam in the machine. The global 
"pedestal shift" variable Pedestal Shift gets subtracted from the tower pedestal, 
in order to center the pedestal-subtracted tower signal at Pedestal Shift. Finally, 
the pedestal gets rounded to the nearest multiple of four and two last bits are 
removed. The last four bits of the result are used, together with the sign, to fill 
the 5 bit pedestal register PED in the FEE. 

During 2003 data taking run, the pedestal subtraction scheme was not yet 
implemented in FEEs, and the HighTower sensitivity was chosen to be HT6 = 3. 
For the 2005 data, the settings were Pedestal Shift = 24 and HT6 = 2, which 
defined ADC™|g^°™' = ADCP^^^,^"'" = 16 x ADCtower and thus aligned the 
HighTower and PatchSum readings at the center of bin 1 in the absense of tower 
physics signal. 

The LUT arrays are prepared in a way that gives a linear response to the 
patch sum in the range from to 62, to allow diagnosing the broken cables by 
observing "all ones" bit pattern 63 during the run. With the nominal setting 
of PedestalShift = 24, each tower contributes a pedestal value 24/16 = 1 to the 
sum, so the LUT is constructed from the following three pieces: 

• Zero, if the sum of 16 towers is below 16 
LUT(s) =0, < s < 16 

• Linear rise in response to the sum in the 6 bit range, excluding 63 
LUT(s) = s - 15, 16 < s < 16 + 63 

• Saturation at 62 

LUT(s) =62, s > 16 + 63 

If a tower is masked out of the PatchSum trigger, the LUT is modified to 
accomodate the loss of its pedestal, in this case it starts rising one count earlier 
than the nominal 16. Therefore, the so-called LUT pedestal (PEDlut) is equal 
to the number of masked towers in a patch. 

The SMD electronics (FEE) board is mounted on the |?7| — 1 side of each 
module. At the FEE board, the amplified cathode strip signals are buffered in 
a switched capacitor array (SCA) before being delivered to external digitizer 
boards outside of the STAR magnet. 
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The signals from the pads of the SMD are amphfied and stored in an analog 
pipeline, composed of switched capacitor arrays, to await the level-0 trigger. 
Upon level-0 trigger, the SMD analog signals are queued with multiplexing 
ratio of 80 : 1 to the 10-bit SMD digitizers. SMD digitized signals are first 
available in STAR level-2 trigger processors in 200 /iS, still well ahead of digital 
information from the TPC. 

The digitizing electronic boards and crates for preshower detector are iden- 
tical to the ones used in the SMD. 
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Summary 



This thesis presents a measurement of neutral pion and eta meson production 

in p + p and d + Au collisions at a ccnter-of-mass energy of y^T^N = 200 GeV, 
measured with STAR detector at Relativistic Heavy Ion Collider (RHIC) in 
Brookhaven National Laboratory (BNL, USA). 

The present neutral pion spectrum complements that of the charged pions 
measured in STAR in the transverse momentum range 0.35 < pt < 10 GeV/c 
and extends up to pr = 17GeV/c. There is a good agreement between the 
neutral and charged pion cross sections in STAR, in spite of very different 
methods and detectors used in the analysis. The neutral pion cross section 
also agrees well with the measurements of PHENIX, another large detector at 
RHIC, and with the theoretical NLO pQCD calculations. 

This thesis also presents the first measurements by STAR of eta meson 
production, which are in agreement with the PHENIX measurements and with 
the mr-scaling assumption. 

Possible medium-induced modifications of particle production in a nucleus- 
nucleus collision, compared with an incoherent superposition of nucleon-nucleon 
collisions, can be observed by measuring the so-called nuclear modification 
factor R. This thesis presents the measurements of RdA, where the neutral 
pion production in d + Au collisions is compared to that in the p + p, and Rcp, 
the comparison between central and peripheral d + Au collisions. Both results 
are in a good agreement with the charged pion measurements previousely done 
by STAR. 

The p + p and d + Au results presented here provide a baseline measurement 
for the future Au 4- Au measurements. These measurements are interesting to 
shed light on quark number scaling in particle production at intermediate pt 
and to study the origin of suppression phenomena at large pt- 
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Samenvatting 



Dit proefschrift beschrijft metingen van neutrale pion en eta meson productie 

in p+p en d + Au botsingen bij een zwaartepuntsenergie van y^s^N = 200 GeV, 
gemeten met de STAR detector bij de Relativistic Heavy Ion Collider (RHIC) 
in Brookhaven National Laboratory (BNL, USA). 

Het neutrale pion spectrum is complementair aan dat van de geladen pioncn 
gemeten met STAR in het transversale impuls gebied 0.35 < pr < 10 GeV/c 
en loopt door tot pt = 17GeV/c. Er is een goede overeenstemmingtussen de 
neutrale en geladen werkzame doorsnede in STAR, ondanks de zeer verschillende 
methoden en detectoren die gebruikt zijn in de analyse. De gemeten werkzame 
doorsnede van de 7r° stemt ook goed overeen met de metingen van PHENIX, 
een andere detector bij RHIC, en met NLO pQCD berekeningen. 

Dit proefschrift beschrijft ook de eerste metingen met STAR van eta meson 
productie, die in overeenstemming zijn met de PHENIX metingen en met de 
mT-schaling aanname. 

De door het medium geinduceerde verandering van deeltjes productie in 
een kern-kern botsing, vergeleken met incoherente superpositie van individuele 
nucleon-nucleon botsingen, kan waargenomen worden door de zogenaamde 
"nuclear modification factor" R te meten. Dit proefschrift beschrijft de metin- 
gen van RdA, waar productie in Au botsingen vergeleken worden met die 
in p + p botsingen, en RcP: dc vcrgelijking tussen centrale en perifere d + Au 
botsingen. Beide resultaten zijn in goede overeenstemming met eerdere metin- 
gen aan geladen pionen met STAR. 

De p+p en Au resultaten die hicr gepresenteerd worden zijn ccn rcfcrcntic 
voor toekomstige Au + Au metingen. Deze metingen zijn intcrcssant omdat ze 
licht werpen op de schaling met het aantal quarks in deeltjes productie bij inter- 
mcdiairc px en voor de studie naar de oorsprong van onderdrukkingsfenomenen 
bij hoge Pt- 
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